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19.  Abstract  -  ( Continued) 


An  interim  report  is  included  ora  study  of  about  100  local  and  regional  events  detected  by  NORESSy 
This  study  is  a  systematic  effort  to  test  analytical  methods  being  developed^at  the  Centei*  for  extracting 
features  from  seismic  waveforms  that  can  be  used  to  identify  regional  phases,  and  for  recognizing  repeated 
events  from  a  single  source  such  as  quarry  blasts.  The  methods  under  evaluation  include  the  extraction  of 
frequency-domain  spectral  parameters,  time-domain  spectral  paramenters,  and  particle  motion  informa¬ 
tion.  This  work  has  particular  application  to  analysis  of  three-component  data  at  GSE  "National  Data 
Centers."  However,  the  work  on  recognizing  repeated  sources  is  also  applicable  to  array  stations,  and  the 
phase  identification  work  is  also  applicable  to  identifying  phases  not  separable  by  conventional  f-k  analysis 
from  array  data  (e.g.,  Sn  vs.  Lg  ). 

Substantial  effort  during  this  quarter"  was  devoted  to  developing  concepts  and  preparing  for  experi¬ 
ments  in  international  exchange  of  seismic  waveform  data.  This  work  included  developing  the  technical 
concepts  for  a  potential  global  seismic  monitoring  system, {with  much  interaction  with  DARPA)  and 
preparing  the  report  (GSE/US/44)  that  was  subsequently  tabled  in  Geneva  at  the  GSE  meeting  in  March. 
Other  GSE  work  included  installing  and  debugging  the  "automatic  association"  program  on  the  Center’s 
GSE  Sun  computer  (in  lieu  of  a  VAX),  and  providing  the  program  to  Australia  for  the  data  center  being 
developed  there  for  GSE  experiments.  Waveform  exchange  experiments,  initiated  during  1986  with  the 
exchange  of  a  few  waveforms,  became  an  important  activity  during  this  quarter  as  described  in  this  report. 

Our  investigation  of  the  sensitivity  of  mi  to  changes  in  t*  have  now  been  extended  to  explosions  at 
twelve  sites.  This  work  continues  to  show  that  mi  is  least  affected  by  At*  if  Q  is  frequency  dependent,  and 
that  the  mh  correction  for  At*  is  a  function  of  magnitude  itself  (as  expected  from  theory).  An  apparent 
dependence  of  magnitude  anomalies  on  explosion  magnitude  has  been  noted  for  Soviet  tests  recorded  at 
NORESS.  Capabilities  of  three-component  stations  to  determine  teleseismic  source-direction  by  polariza¬ 
tion  analysis  are  nearing  completion,  and  show  that  a  useful  analytical  tool  for  GSE  National  Data  Centers 
is  available. 

The  Center  has  been  heavily  involved  in  an  examination  of  its  database  management  system  and  the 
development  of  extensions  to  handle  new  types  of  data— principally  data  from  seismic  arrays  and  data  anti¬ 
cipated  from  forthcoming  GSE  experiments.  Working  jointly  with  Science  Horizons,  Inc.,  and  S AIC’s  La 
Jolla-based  geophysics  group,  the  Center  plans  to  propose  a  revision  of  its  version  2.7  database  structure 
during  the  next  quarter.  Major  work  on  building  databases  containing  NORESS  data  has  continued,  and 
progress  reports  on  this  and  other  database  development  are  included. 
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1.  FOREWORD 


Work  continued  during  this  reporting  period  on  the  assessment  of  capabilities  of 
seismic  networks  and  the  design  of  networks  to  achieve  specified  performance  characteris¬ 
tics.  Results  reported  in  Section  2  show  that  a  network’s  location  accuracy  approaches  a 
near-maximum  when  stations  are  distributed  over  an  azimuthal  sector  of  about  180 
degrees  and  that  good  depth  determination  requires  stations  within  about  15  degrees  of 
the  focus  of  the  event.  The  latter  suggests  that  global  networks  should  include  100  or 
more  stations  if  accurate  depth  estimation  is  important.  As  in  previous  reports,  this  work 
continues  to  be  directed  toward  networks  of  the  type  being  considered  by  the  Group  of 
Scientific  Experts  (GSE). 

Section  3  presents  an  interim  report  on  a  study  of  about  100  local  and  regional 
events  detected  by  NORESS.  This  study  is  a  systematic  effort  to  test  analytical  methods 
being  developed  at  the  Center  for  extracting  features  from  seismic  waveforms  that  can  be 
used  to  identify  regional  phases,  and  for  recognizing  repeated  events  from  a  single  source 
such  as  quarry  blasts.  The  methods  under  evaluation  include  the  extraction  of  frequency- 
domain  spectral  parameters,  time-domain  spectral  paramenters,  and  particle  motion 
information.  This  work  has  particular  application  to  analysis  of  three-component  data  at 
GSE  "National  Data  Centers."  However,  the  work  on  recognizing  repeated  sources  is  also 
applicable  to  array  stations,  and  the  phase  identification  work  is  also  applicable  to  identi¬ 
fying  phases  not  separable  by  conventional  f-k  analysis  from  array  data  (e.g.,  Sn  vs.  Lg  ). 

Substantial  effort  during  this  quarter  was  devoted  to  developing  concepts  and 
preparing  for  experiments  in  international  exchange  of  seismic  waveform  data.  This  work 
included  developing  the  technical  concepts  for  a  potential  global  seismic  monitoring  sys¬ 
tem  (with  much  interaction  with  DARPA)  and  preparing  the  report  (GSE/US/44)  that 
was  subsequently  tabled  in  Geneva  at  the  GSE  meeting  in  March.  Other  GSE  work 
included  installing  and  debugging  the  "automatic  association"  program  on  the  Center’s 
GSE  Sun  computer  (in  lieu  of  a  VAX),  and  providing  the  program  to  Australia  for  the 
data  center  being  developed  there  for  GSE  experiments.  Waveform  exchange  experi¬ 
ments,  initiated  during  1986  with  the  exchange  of  a  few  waveforms,  became  an  important 
activity  during  this  quarter;  details  are  given  in  Section  4  of  this  report. 

Progress  reports  on  a  number  of  other  activities  at  the  Center  are  given  in  Section  5. 
Our  investigation  of  the  sensitivity  of  to  changes  in  t*  have  now  been  extended  to 
explosions  at  twelve  sites.  This  work  continues  to  show  that  m ^  is  least  affected  by  At* 
if  Q  is  frequency  dependent,  and  that  the  correction  for  At*  is  a  function  of  magni¬ 
tude  itself  (as  expected  from  theory).  An  apparent  dependence  of  magnitude  anomalies 
on  explosion  magnitude  has  been  noted  for  Soviet  tests  recorded  at  NORESS.  Capabili¬ 
ties  of  three-component  stations  to  determine  teleseismic  source-direction  by  polarization 
analysis  are  nearing  completion,  and  show  that  a  useful  analytical  tool  for  GSE  National 
Data  Centers  is  available. 


t 


The  Center  has  been  heavily  involved  in  an  examination  of  its  database  management 
system  and  the  development  of  extensions  to  handle  new  types  of  data-principally  data 
from  seismic  arrays  and  data  anticipated  from  forthcoming  GSE  experiments.  Working 
jointly  with  Science  Horizons,  Inc.,  and  S AIC’s  La  Jolla-based  geophysics  group,  the 
Center  plans  to  propose  a  revision  of  its  version  2.7  database  structure  during  the  next 
quarter.  Major  work  on  building  databases  containing  NORESS  data  has  continued,  and 
progress  reports  on  this  and  other  database  development  is  also  contained  in  Section  5. 
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2. 


NETWORK  CAPABILITY  AND  DESIGN 


2.1.  A  NOTE  ON  NETWORK  DESIGN:  STATION  DISTRIBUTION 
AND  EPICENTER  ERROR 


2.1.1.  Introduction 

Determination  of  the  epicenter  of  a  detected  seismic  event  is  one  important  function 
for  seismological  test  ban  monitoring.  For  example,  a  precise  epicenter  determination  is 
important  for  an  event  the  nature  of  which  can  not  readily  be  established.  This  is  the 
case  particularly  if  it  occurs  close  to  territorial  boundaries  or  if  such  an  event  is  selected 
for  on-site  inspection. 

The  quality  of  an  epicenter  estimate  depends  on  several  factors  one  of  which  is  the 
geometry  of  the  detecting  stations  in  relation  to  the  epicenter.  A  symmetrical  distribu¬ 
tion  of  the  stations  around  the  azimuth  from  the  event  is  often  considered  desirable.  This 
is,  however,  often  not  possible  to  achieve  in  practice.  In  this  note  we  attempt  to  describe 
how  uncertainties  in  epicenter  determination  vary  with  the  relative  location  of  the  detect¬ 
ing  seismic  stations.  The  analysis  is  illustrated  with  calculated  location  errors  for 
hypothetical  station  networks. 

2.1.2.  Model 

Estimation  of  the  epicenter  (latitude  and  longitude)  of  an  event  is  part  of  the  seismic 
event  location  procedure,  which  also  simultaneously  gives  estimates  of  origin  time  and 
focal  depth.  In  this  note  the  analysis  is,  however,  limited  to  the  joint  marginal  confidence 
region  for  the  estimated  epicenter  (Flinn,  1965). 

We  introduce  a  Cartesian  coordinate  system,  X=(x1,x2),  with  origin  at  the  true 
epicenter.  The  epicenter  determination  is  assumed  to  be  based  on  first-arrival  time  read¬ 
ings  at  n  stations,  which  are  located  at  distances  Ak  and  azimuths  4>k  from  the  epicenter. 
The  arrival  time  residuals  or  data  errors  are  assumed  to  be  normally  distributed  and  are 
characterized  by  an  a  priori  standard  deviation,  a.  These  data  errors  include  reading 
errors  at  the  stations  as  well  as  uncertainties  of  the  travel  time  model.  Since  standard 
error  is  usually  dominated  by  uncertainties  in  the  travel  time  model,  it  is  assumed  here 
that  all  observations  have  the  same  a. 

The  joint  marginal  confidence  region  for  the  estimated  epicenter  is  given  by 

a-‘-XTSXix'(2) 
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The  elements  of '  ie  symmetrical  covariance  matrix  S  are  determined  by  the  relative  loca¬ 
tions  of  the  n  detecting  stations  as  described  by  their  azimuths  from  the  source 

dT 

(  ^t;&=l,2,...,n  )  and  their  slowness  values  ( - ):  . 

dAk 

dT  2  , 

5n=E( - )  ■cosfyj 

dAk 
dT  2 

522=E( - ) 

dAk 
dT  2 

512=E( - )  -cos^J-sin^) 

dAk 

In  order  to  simplify  the  expressions  above  for  the  elements  of  the  matrix  S  is  assumed 
that  the  epicenter  of  the  event  is  located  at  the  north  pole.  The  confidence  region  is  an 
ellipse  and  the  orientation  and  lengths  of  the  axes  (  a  and  b)  are  thus  determined  by  the 
elements  of  the  matrix  S. 

We  will  use  the  area  (i.e,,  ira'b)  and  the  square  root  of  the  area  of  this  confidence 
ellipse  as  measures  of  the  location  error.  Since  we  consider  location  errors  in  contexts  like 
on-site  inspection,  the  area  is  used  rather  than  the  major  axis  as  a  measure  of  the  location 
error. 

It  can  be  shown  that  the  area  of  the  confidence  ellipse  is  inversely  proportional  to 
the  determinant  of  the  matrix  S: 

*-a-b=*x2{2y<T2/{Sn-S22-Sl2) 

This  formula  is  180  degrees  symmetric  with  regard  to  the  azimuthal  station  distribution. 
In  other  words  if  a  station  is  moved  180  degrees  around  the  source  azimuth  the  location 
error  will  still  be  the  same.  This  means  that  a  360  degree  coverage  of  stations  around  the 
source  azimuth  will  not  give  a  smaller  location  error  than  a  coverage  of  180  degrees.  This 
conclusion  is  based  on  the  particular  assumptions  about  the  model  used  here;  in  practice, 
due  to  lateral  variation  of  the  velocity  in  the  source  region  systematic  azimuthal  bias  may 
be  introduced  with  stations  in  only  a  180  degree  sector  around  the  source. 

If  we  assume  that  all  stations  are  at  the  same  distance  from  the  source  with  equal 
slowness  values  (i.e.,  dT/dAk=P),  the  determinant  of  S  can  be  written  as: 

2 

I  S|  =p{  Ecos2(^*)|I  Esin2(^*)l  4  Esin(^*)'cos(^*)| 


According  to  Cauchy’s  inequality  the  determinant  is  always  greater  than  zero  except 
when:  sin(<^)  =  c  •cos(<^i),  for  all  k,  where  e  is  a  constant. 

In  order  to  study  distributions  with  a  maximum  value  of  the  determinant  we  write 
xt=sin(^fc)  and  this  gives  the  following  expression  for  the  determinant: 

|  Xftjj  £(i-**)}  4 

Rather  than  attempting  to  find  the  xk  values  for  a  given  n  that  maximize  the  deter¬ 
minant,  we  assume  that  the  xk  are  realizations  of  a  stochastic  variable  £  with  the  range  of 
variation  (-1,1).  We  also  replace  the  summations  with  expectation  operators,  E,  and  get: 

We  are  interested  in  finding  the  distribution  of  £  that  maximizes  the  expression  above. 
For  simplicity  we  assume  that  £  is  uniformly  distributed  over  the  interval  a,  (3.  In  this 
case  the  determinant  can  then  be  written  in  a  closed  form: 


For  a  and  (3  we  have: 


Figure  1  shows  the  relative  location  error  as  a  function  of  azimuthal  coverage  for  a  "uni¬ 
form”  distribution  as  defined  above.  The  curve  suggests  that  error  decreases  rapidly  with 
increasing  coverage  at  small  values  of  the  coverage  and  a  coverage  of  90  degrees  gives  an 
error  about  20  percent  larger  than  that  with  a  coverage  of  180  degrees. 

2.1.3.  Some  Examples  of  Uniform  Distributions 

In  order  to  study  the  effect  of  distance  and  azimuth  coverage  on  location  errors  we 
calculate  the  location  error  for  some  hypothetical  networks,  which  are  distributed  along 
circular  arcs  around  the  epicenter  at  various  epicentral  distances,  i.e., 

dT 


for  all  stations,  k. 


0.00 


30.00  60.00  90.00  120.00  150.00 

Azimuth  Coverage  (degrees] 


180. 00 


Figure  1.  The  curve  shows  the  relative  location  error  as  a  function  of  azimuth  for  a  uni 
form  station  distribution  represented  by  a  stochastic  variable  as  described  in  the  text. 


We  assume  that  the  event  occurs  at  a  depth  of  1  km,  and  calculate  errors  for  networks  on 
arcs  with  a  radius  of  10,  20,  30,  90  degrees.  The  length  of  the  arc  or  the  azimuth 

covered  vary  from  22.5  to  270  degrees.  The  number  of  stations  in  the  networks  are 
n=4,8,12,16,...,32. 

In  this  case  the  area  of  the  error  ellipse  simplifies  to  the  formula: 

*  S|  /nP(A) 

here  the  elements  of  the  matrix  S  do  not  contain  the  slowness  values.  We  note  that  the 
square  root  of  the  area  of  the  confidence  ellipse  is  inversely  proportional  to  the  slowness 
value  and  directly  proportional  to  the  a  priori  standard  deviation. 

In  other  words,  the  area  of  the  location  error  is  direct  proportional  to  the  number  of 
stations,  n. 

Examples  of  some  hypothetical  networks  with  associated  error  ellipses  are  shown  in 
Figure  2.  The  figure  includes  six  networks  distributed  on  circular  arcs  at  10,  20,  30,  50, 
70,  and  90  degrees  with  4,  8,  12,  16,  20,  and  24  stations  with  azimuthal  coverage  of  22.5, 
45.0,  90.0,  135.0,  180.0,  and  270.0  degrees.  For  example  the  network  at  50  degrees  con¬ 
sists  of  16  stations  covering  an  azimuthal  sector  of  135  degrees.  The  associated  error 
ellipses  are  to  the  right  in  the  figure.  The  two  largest  ellipses  corresponds  to  the  smallest 
networks  with  4  and  8  stations.  In  can  be  seen  that  there  is  no  dramatic  difference  in  the 
area  of  the  other  four  ellipses,  which  are  associated  with  networks  covering  an  azimuthal 
sector  of  90  degrees  or  more. 

The  effect  of  azimuthal  coverage  on  the  location  error  can  also  be  seen  in  Figure  S, 
which  shows  the  error  ellipses  for  4  and  32  station  networks  both  with  the  stations  at  60 
degrees  epicentral  distances  and  at  varying  azimuthal  coverage  from  22.5  to  270.0  degrees. 
In  this  case  the  effect  of  distance  is  kept  fixed  and  the  effect  of  number  of  stations  is  con¬ 
trolled  by  looking  at  extreme  values.  Again,  it  can  be  seen  that  there  is  no  dramatic 
difference  between  the  error  ellipses  for  networks  with  a  90  degree  or  more  coverage  for  a 
given  number  of  stations.  As  expected,  the  number  of  stations  affect  the  overall  size  of 
the  error  ellipses. 

Figure  4  shows  error  ellipses  for  networks  with  a  given  azimuthal  coverage  (90 
degrees)  with  the  stations  (in  all  12)  at  different  distances  from  10  to  90  degrees  in  the 
left  diagram,  and  with  different  numbers  of  stations  (4  to  32)  at  a  given  distance  (60 
degrees)  in  the  right  diagram.  The  areas  of  these  ellipses  vary  with  the  slowness  as  a 
function  of  epicentral  distance  and  with  the  square  root  of  the  number  of  stations. 

The  location  error  as  a  function  of  azimuthal  coverage,  epicentral  distance,  and 
number  of  stations  for  concentric  networks  are  outlined  by  Figures  5-  7.  In  Figure  5  the 
location  error  is  shown  for  a  12  station  network  with  azimuthal  coverages  from  22.5 
degrees  and  distances  10  to  90  degrees.  The  location  error  drops  about  a  factor  of  two 
as  the  azimuthal  coverage  increases  from  22.5  to  120  or  more  degrees.  The  curves  in 


•Vi 

& 


m 

[t'j 


•N 

,1  v. 


\«d 

v* 


» 


(  a 

v2 


** 


n 


Figure  2.  Examples  of  concentric  station  networks  with  4,  8,  12,  16,  20,  and  24  stations 
at  epicentral  distances  10,  20,  30,  50,  70,  and  90  degrees,  respectively,  and  with  azimuthal 
sector  coverages  of  22.5,  45.0,  90,  135.0,  180,  and  270.0  degrees  (to  the  left).  The  figure 
also  shows  the  associated  error  ellipses  of  the  epicenter  estimates  (to  the  right).  The  two 
largest  ellipses  are  associated  with  the  networks  of  4  and  8  stations. 


Figure  8.  The  figure  shows  the  error  ellipses  for  concentr&l  networks  at  60  degrees  epi- 
central  distance  and  with  4  (to  the  left)  and  32  (to  the  right)  stations  covering  azimuthal 
sectors  of  22.5,  45.0,  90.0,  135.0,  180.0,  and  270.0  degrees.  The  two  largest  ellipses  in 
each  diagram  corresponds  to  the  networks  with  22.5  and  45.0  azimuthal  degree  coverage. 
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Figure  4.  Error  ellipses  of  concentral  station  networks  with  90  degrees  azimuthal  cover¬ 
age  and  with  12  stations  at  varying  epicentral  distances  (10-90  degrees,  to  the  left)  and 
with  varying  number  of  stations  (4-32)  at  60  degree  epicentral  distance  (to  the  right). 
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Figure  5.  Location  error  as  a  function  of  azimuthal  coverage  for  concentric  networks 
with  12  stations  at  epicentrai  distances  10  to  90  degrees. 
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Figure  5  also  shows  an  even  more  significant  decrease  (about  a  factor  of  three)  in  error  as 
the  epicentral  distance  of  the  concentric  network  decreases  from  90  to  10  degrees. 

The  distance  effect  is  also  shown  in  Figure  6,  whereas  the  error  is  shown  for  a  net¬ 
work  of  12  stations  at  azimuthal  coverages  from  22.5  to  180.0. 

Finally,  in  Figure  7  the  error  is  shown  as  a  function  of  number  of  stations  for  concen¬ 
tric  networks  with  90  degrees  azimuthal  coverage  at  distances  from  10  to  90  degrees.  A 
rather  small  reduction  in  location  error  is  obtained  for  15  stations  or  more. 

2.1.4.  Concluding  Remarks 

In  this  note  we  analyze  epicenter  errors  for  hypothetical  network  with  regard  to 
number  of  stations,  azimuthal  distribution,  and  epicentral  distances  with  a  standard  sta¬ 
tistical  model.  We  use  the  area  of  the  joint  marginal  confidence  ellipse  as  a  measure  of 
the  location  error. 

The  errors  due  to  azimuth  distribution  are  periodic  with  a  period  of  180  degrees,  i.e., 
an  azimuthal  coverage  of  360  degrees  around  the  epicenter  does  not  provide  smaller  errors 
than  a  coverage  of  180  degrees.  For  networks  with  a  uniform  distribution  along  an  epi¬ 
central  arc  around  the  event  epicenter  a  coverage  of  90  degree  gives  errors  which  are  not 
dramatically  larger  than  those  of  networks  with  a  complete  180  degree  coverage.  This 
conclusion  is  based  on  the  particular  assumptions  about  the  model  used;  in  practice, 
lateral  velocity  variations  in  the  source  region  may  introduce  systematic  azimuthal  bias 
with  stations  covering  a  sector  of  180  degrees  or  less. 

The  epicentral  error  is  inversely  proportional  to  the  slowness  for  networks  with  sta¬ 
tions  at  equal  epicentral  distance  from  the  event. 

The  area  of  the  location  error  is  for  such  networks  inversely  proportional  to  the 
number  of  stations  of  the  network. 
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2.2.  A  NOTE  ON  NETWORK  DESIGN:  STATION  DISTRIBUTION 
FOR  MINIMUM  VARIANCE  OF  DEPTH  ESTIMATES 


2.2.1.  Introduction 

Determination  of  focal  depth  is  one  important  function  of  a  seismological  station 
network  for  test  ban  monitoring.  Accurate  depth  estimates  can  be  obtained  on  the  basis 
of  surface  reflected  phases,  pP  and  aP.  Unfortunately  the  frequency  of  reporting  of  sur¬ 
face  reflections  is  usually  quite  low,  and  depth  phases  are  also  difficult  to  identify,  in  par¬ 
ticular  for  shallow  events.  Even  if  the  prospects  for  comprehensive  use  of  surface  reflec¬ 
tions  for  depth  determination  appear  rather  limited,  there  are  also  studies  that  indicated 
that  significantly  improved  observability  may  be  obtained  if  waveform  data  from  a  net¬ 
work  of  stations  can  be  analyzed  simultaneously  (Roy,  1984).  Only  systematic  future 
attempts  of  utilising  such  a  network  approach  with  waveform  data  for  detection  of  surface 
reflections  can  confirm  these  promising  tentative  results.  For  the  time  being  it  appears 
that  first  arrival  times  of  P- waves  will  be  the  most  important  information  available  for 
depth  determinations,  in  particular  for  small  magnitude  events  for  which  waveform  data 
is  nothing  more  than  a  small  amplitude  barely  above  the  background  noise.  In  this  note 
we  discuss  the  relative  location  of  stations  and  seismic  events  for  accurate  depth  estima¬ 
tion  based  on  first  arrival  times. 

2.2.2.  Model 
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We  use  a  simplified  approach  to  study  the  accuracy  of  depth  estimates  as  a  function 
of  geographical  distribution  of  stations  (Lomnitz,  1977).  It  is  assumed  that  the  epicenter 
of  the  event  is  known,  the  focal  depth  will  be  estimated  and  that  n  stations  have  reported 
arrival  times,  f  .,  which  all  have  an  error  with  the  same  standard  deviation  <r.  Forming 
the  residual,  tei,  between  the  observed  and  the  calculated,  (*„■),  arrival  times,  the 

following  relations  are  obtained  between  the  residuals,  the  adjustments  in  depth  ,  dh,  and 
in  origin  time  dr. 


r{  =  dr+dh 
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This  can  be  regarded  as  a  straightforward  linear  regression  relation  and  the  variance  of 
the  estimated  dh  can  be  written  as: 


V{dh)= 
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For  standard  travel  time  models  the  derivative - is  a  function  of  epicentral  distance  A 

dh 

only  for  a  given  depth  h.  Examples  of  this  derivative  as  a  function  of  distance  for  the 
Jeffreys-Bullen  travel  time  model  for  focal  depths  at  10,  30,  50,  100,  200  km  are  shown  in 
Figure  1.  The  estimated  curves  were  obtained  by  a  bicubic  two-dimensional  spline 
approximation  to  Jeffreys-Bullen  travel  time  tables.  The  curves  for  depths  at  30  km  and 
more  are  monotonically  decreasing  with  distance  with  most  of  the  decrease  occurring  at 

distances  less  than  20  degrees.  Since  the  range  of  variation  of  -  is  limited,  the 

dh 

minimum  variance  of  dh  is  obtained  if  all  the  stations  are  located  at  the  boundaries  of  the 

distance  interval  over  which  they  are  distributed.  In  other  words,  if  there  are  n  stations 

(assume  even  number)  that  can  be  located  anywhere  within  the  distances  interval  Amin  to 

A _ _  the  minimum  variance  is  obtained  if  half  of  the  stations  are  located  at  A  and 

max’  nun 

the  other  half  at  A _ ,  and  the  standard  deviation  of  dh  becomes: 

max7 


<r(dh)  =  2  -<r 


x^idt/dh(A.)-dt/dh(A)) 


2.2.3.  Examples 

The  formula  above  for  the  variance  of  the  depth  estimate  with  stations  symmetri¬ 
cally  placed  at  the  ends  of  the  range  of  intervals  covered  has  been  used  for  a  network  of 
four  stations,  two  placed  at  either  end  of  the  covered  distance  interval  in  Figure  2.  Two 
of  the  stations  are  supposed  to  be  at  80  degrees  distance  and  the  standard  deviation  is 
shown  as  a  function  of  the  distance  of  the  two  other  stations  (minimum  distance)  for 
events  at  30  and  200  km  depths.  The  standard  deviation  starts  to  increase  significantly 
as  the  minimum  distance  becomes  larger  than  15  degrees  or  so.  In  fact,  it  is  increased 
four  times  as  it  goes  from  10  to  30  degrees.  This  illustrates  the  importance  of  one  or  more 
detecting  stations  at  local  or  regional  distance  in  a  network  based  primarily  on  teleseismic 
detection. 

If  we  assume  that  a  network  has  a  fairly  equal  station  spacing  on  the  globe  we  can 
ask  for  any  given  location  what  is  the  probability  that  stations  within  say  15  degrees  do 
not  detect  an  event  and  four  or  more  stations  at  distances  beyond  15  degrees  do  detect 
the  event.  This  probability  is  denoted  P{aiargt)  below.  In  this  case  the  standard  devia¬ 
tion  of  the  estimated  depth  would  be  fairly  large,  and  a  low  probability  is  desirable.  If  we 
assume  a  network  of  N  stations,  the  number  of  stations  within  A  degrees  is: 
n(A)  =  (l— cos(A))/2,iV  for  equal  station  spacing.  The  detection  probability  at  station  i 
as  a  function  of  magnitude  is  usually  described  by  a  normal  distribution  function  and  by 


the  station  threshold  T-  and  associated  standard  deviation  a.:  p(m)  =  $j  (m— T( 
The  probability  mentioned  above  can  be  written  as: 

pK„.)=p<^4)p(’>=°) 
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Figure  1.  The  travel  time  depth  derivative,  dt/dh ,  as  a  function  of  epicentral  distance 
derived  from  the  Jeffrey-Bullen  P-travel  time  tables. 
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Figure  2.  Standard  deviation  of  estimated  focal  depth  for  events  detected  by  four  sta¬ 
tions.  Two  stations  are  80  degrees  from  the  epicenter,  and  the  other  two  are  at  various 
distances  as  indicated  by  the  horizontal  axis  of  the  figure.  Notice  that  the  standard  devia¬ 
tion  increases  rapidly  as  this  minimum  distance  increases  from  10  to  25  degrees. 
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where  £  and  rj  are  binomial  stochastic  variables  for  the  number  of  detecting  stations 
between  15-80  and  less  than  15  degrees,  respectively.  Assuming  the  thresholds  TV  are 
4.72  and  5.0  for  station  distance  less  than  15  and  between  15-80  degrees,  respectively,  and 
<7=0.40  we  get  the  probability  curves  of  Figure  S  as  a  function  of  magnitude,  m,  for  net¬ 
works  with  N= 50  to  1000  stations.  The  recently  suggested  GSE  concept  network 
includes  50  fairly  uniformly  distributed  stations.  The  curves  suggest  that  a  significant 
reduction  of  P(^iarft)  is  only  obtained  with  increasingly  large  number  of  stations,  several 
hundred  or  more  for  a  complete  coverage  of  the  globe.  If  desirable  coverage  is  limited  to 
the  continents  these  numbers  would,  however,  be  reduced  accordingly. 

2.2.4.  Concluding  Remarks 

In  this  note  we  study  the  standard  deviation  of  estimated  focal  depth  on  the  basis  of 
first-arrival  times. 

It  is  shown  this  standard  deviation  is  entirely  determined  by  the  distance  range 
covered  but  detecting  stations.  Coverage  limited  to  the  teleseismic  range  30-80  degrees 
will  give  comparatively  large  standard  deviations,  whereas  a  significant  improvement  is 
obtained  if  observations  at  15  degrees  or  closer  can  be  utilized. 

An  example  with  idealized  hypothetical  networks  suggests  that  global  station  net¬ 
works  need  to  consist  of  a  large  number  of  stations,  several  hundred,  in  order  to  achieve  a 
comparatively  low  standard  deviation  of  depth  estimates  for  events  throughout  the  world. 
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2.3.  A  NOTE  ON  NETWORK  DESIGN:  SHIELD  AND  PLATFORM 
NETWORKS 


Seismological  stations  located  on  shields  and  platforms  are  often  observed  to  have 
lower  event  detection  thresholds  than  those  for  stations  located  in  tectonically  active 
regions.  The  simplicity  of  shield  and  platform  structures  also  frequently  results  in  com¬ 
paratively  simple  recorded  seismic  signals  which  facilitates,  for  example,  identification  of 
depth  phases  and  waveform  analysis,  which  has  become  an  important  new  element  of  the 
GSE-system.  With  these  advantages  to  be  accrued  from  shield  and  platform  stations,  we 
are  exploring  the  capabilities  of  global  networks  composed  of  stations  primarily  located  on 
shields  and  platforms.  This  note  summarized  the  status  of  this  investigation. 

Shield  and  platform  areas  that  have  had  seismological  stations  located  in  them  are 
marked  in  Figure  1,  which  is  based  on  a  station  list  including  about  4000  stations  com¬ 
piled  at  the  Center  for  Seismic  Studies,  which  in  turn  is  primarily  based  on  compilations 
by  the  USGS  (USGS,  1985).  There  are  58  and  75  segments  (five  by  5  degrees  in  latitude 
and  longitude  according  to  regionalization  by  Jordan,  1984)  with  one  or  more  stations  for 
shield  and  platform  regions,  respectively.  In  all,  this  gives  133  segments  with  one  or  more 
stations.  This  number  can  be  compared  with  the  number  (20)  of  shield  and  platform  sta¬ 
tions  that  participated  in  the  GSETT. 

Calculations  of  event  detection  capabilities  and  other  seismological  monitoring  func¬ 
tions  were  made  for  several  hypothetical  networks  including  one  consisting  of  133  stations 
located  at  the  segments  where  seismic  stations  have  been  in  operation.  For  the  sake  of 
comparison,  calculations  were  made  for  networks  with  25  -  175  stations  more  uniformly 
distributed  over  the  shield  and  platform  areas,  and  for  a  hypothetical  version  of  the 
GSETT  network  (total  of  69  stations)  and  the  GSE  concept  network  (GSE/US/44,  1987). 
The  latter  consists  of  50  stations,  eight  of  which  have  arrays.  The  same  station  parame¬ 
ters  were  used  for  all  networks  (the  only  exceptions  were  the  eight  array  stations  of  the 
concept  network,  which  were  given  a  short-period  noise  amplitude  of  2  nm  instead  of  10 
nm).  Notice  that  the  noise  amplitudes  assumed  for  the  GSE  concept  network  here  are  a 
factor  of  two  larger  than  those  used  in  earlier  calculations  (GSE/US/44,  1987).  The  cal¬ 
culations  assumed  only  teleseismic  observations,  i.e.,  use  of  regional  phases,  which  are 
particularly  important  for  depth  estimation  was  not  included. 

Examples  of  calculated  event  detection  thresholds  are  shown  in  Figures  2  and  S. 
The  contours  of  the  90%  mh-  threshold  for  the  GSE  concept  network  are  compared  with  a 
99-station  network  on  shields  and  platforms  in  Figure  2.  The  threshold  for  the  concept 
network  varies  much  less  (4. 4-4. 6)  than  that  for  the  99-station  network  (4. 2-4. 8).  How¬ 
ever,  the  threshold  of  the  99-station  network  is  below  4.4  for  most  of  the  continental 
areas.  Figure  S  summarizes  the  thresholds  for  the  networks.  The  figure  suggests  that 
no  significant  lowering  of  the  world  average  threshold  is  achieved  by  increasing  the 
number  of  stations  beyond  75  or  so  for  shield/platform  networks.  The  concept  network 
has  a  fairly  small  range  of  variation,  as  would  be  expected  from  Figure  2.  The  GSETT 
network  has  a  larger  range  of  variation,  primarily  because  of  the  concentration  of  stations 
in  Europe.  The  GSETT  network  also  has  lower  maximum  threshold  than  that  of  75 
shield/platform  stations,  owing  to  some  oceanic  stations. 
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Figure  1.  The  map  shows  the  shield  and  platform  segments  on  which  one  or  more 
seismological  stations  are  or  have  been  operating  (indicated  as  filled  rectangles). 
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Figure  2.  Contours  of  the  90%  m6-detec.tion  threshold  for  the  GSE  concept  network 
(top)  and  99-station  shield/platform  network  (bottom).  P-signal  detections  at  four  or 
more  stations  are  used  as  event  detection  criterion. 
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Figure  S.  Detection  threshold,  m^,  as  a  function  of  number  of  stations  for 
shield/platform  networks.  The  diagram  shows  mean,  range  of  variation  (minimum  and 
maximum),  standard  deviations  (plus  and  minus  mean)  for  each  network  as  indicated  for 
the  175-station  network  to  the  right  in  the  diagram.  Results  are  also  included  for  the 
GSE  concept,  the  GSETT  networks  and  a  hypothetical  133-station  network  of  existing 
stations  as  marked  in  the  figure.  The  GSE  concept  network  which  includes  50  stations 
have  been  displayed  slightly  to  the  right  in  the  diagram  for  sake  of  clarity. 
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In  conclusion,  the  shield  and  platform  networks  seem  to  have  nonuniform  capabili¬ 
ties  on  a  global  scale  for  basic  seismological  functions.  This  is  not  as  pronounced  for 
event  detection  and  Af#-determination  and  is  comparable  to  other  existing  networks,  but 
is  very  clear  for  location  and  depth.  There  seem  to  be  reasonable  capabilities  on  the  con¬ 
tinents,  provided  sufficient  number  of  uniformly  distributed  stations  would  be  available. 
It  appears  that  a  network  with  reasonably  uniform  coverage  cannot  be  based  entirely  on 
stations  in  shield  and  platform  areas.  Such  a  network  would,  however,  be  sufficient  for 
uniform  coverage  of  the  territorial  areas  on  the  major  continents.  The  GSE  concept  net¬ 
work,  including  50  stations,  provides  fairly  uniform  global  coverage  for  event  detection, 
location,  and  M^-determination.  The  depth  estimation  capability  is,  however,  more  lim¬ 
ited  in  particular  for  shallow  focus  explosions.  Some  of  the  limitations  of  depth  estimation 
capability  may  be  improved  by  use  of  waveform  data  and  use  of  local  and  regional  phases. 


Hans  Israelsson 
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3.  RESEARCH  TO  IMPROVE  ANALYSIS  OF  REGIONAL  SEISMIC 
DATA 


3.1.  CHARACTERISTICS  OF  REGIONAL  SEISMIC  PHASES 
RECORDED  BY  THE  NORESS  ARRAY 


3.1.1.  Introduction 

In  previous  technical  reports,  we  have  presented  papers  outlining  the  development  of 
a  number  of  processing  methods  for  the  extraction  of  signal  parameters  from  NORESS 
data.  These  analyses  include  frequency-domain  spectral  parameters,  time-domain  autore¬ 
gressive  (AR)  parameters,  and  particle-motion  information.  These  methods  are  now 
being  used  at  the  Center  to  test  their  applicability  to  the  problem  of  regional  wave  propa¬ 
gation  and  seismic  discrimination,  in  particular  the  discrimination  of  small  chemical 
explosions  from  natural  microearthquakes  (Ai^<3.5).  The  dataset  we  are  testing  includes 
approximately  100  events  in  Scandinavia  and  the  western  Soviet  Union.  To  date,  approx¬ 
imately  50  events  have  been  analyzed  in  detail.  In  this  paper,  we  present  some  prelim¬ 
inary  results  of  our  study.  Our  goal  is  to  review  the  characteristics  of  the  regional  seismic 
phases  Pn,  Sn,  and  L  .  We  will  focus  on  the  frequency  content,  Q  estimates,  and  the 
three-component  characteristics  of  these  phases. 

3.1.2.  Signal  Processing 

Power  spectral  estimates  of  the  NORESS  data  are  performed  interactively  with  a 
Sun  Microsystems  workstation.  The  central  (NRAO)  channel  is  displayed  in  a  number  of 
frequency  bands  and  the  times  of  Pn,  Sn,  and  Lf  are  chosen.  Then  these  time  picks  are 
used  to  retrieve  the  data  from  the  25  vertical  channels.  To  obtain  an  initial  estimate  of 
the  peak  frequency  and  shape  of  the  spectra,  the  multichannel  data  are  modelled  as  a 
damped  harmonic  oscillator  (DHO)  by  a  method  described  by  Dysart  (1986).  The  peak 
frequency  and  damping  factor  of  the  equivalent  DHO  are  useful  in  describing  the  spectral 
content  of  a  signal  in  terms  of  the  least  number  of  physically  significant  parameters. 
Power  spectral  estimates  are  also  determined  for  the  whole  array  by  averaging  the  spectra 
from  each  vertical  channel.  This  averaging  procedure  reduces  the  variance  in  spectral 
amplitudes  and  results  in  a  spectral  estimate  similar  to  a  conventional  periodogram.  The 
peak  frequency  is  obtained  from  the  instrument-corrected  velocity  spectrum.  Spectral 
amplitudes  above  a  preset  signal-to-noise  ratio  are  then  fit  to  an  idealized  source  function 
to  obtain  estimates  of  the  attenuation  factor  Q,  and  the  corner  frequency. 

Particle-motion  information  is  extracted  using  all  the  three-component  sensors  in  the 
array  by  the  method  described  by  Jurkevics  (1986a,  1986b).  In  this  approach,  the  covari¬ 
ance  matrices  from  all  three-component  sensors  are  computed  in  the  time-domain  and 
averaged  together,  and  the  usual  eigenproblem  is  solved  to  obtain  the  polarization  ellipse 
as  a  function  of  time.  Particle  motions  are  characterized  by  various  attributes  of  the 
polarization  ellipse.  In  this  analysis  a  time  pick  was  made  interactively  on  each  of  the 
regional  phases  of  a  seismogram  and  several  particle-motion  attributes  extracted  in  a 
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short-time  segment  in  fixed  frequency  bands  centered  at  4,  6  and  10  Hz.  The  time  seg¬ 
ments  vary  in  length  in  order  to  isolate  the  most  stable  particle  motions  characteristic  of 
each  phase. 

3.1.3.  Pn 

For  all  of  the  events  examined,  the  peak  frequency  of  the  Pn  wave  was  greater  than 
or  equal  to  the  peak  frequency  of  either  5n  or  Lf.  Magnitudes  of  the  events  studied  range 
from  1.8  to  4.8  (Af^)  with  PR  peak  frequencies  ranging  from  2.0  to  10.0  Hz.  The  varia¬ 
tion  of  Pn  peak  frequency  versus  magnitude  does  not  appear  to  follow  a  systematic  pat¬ 
tern,  but  Pn  peak  frequencies  and  corner  frequencies  both  show  a  dependence  on  distance. 
Peak  and  corner  frequencies  for  close  events  are  generally  higher  than  those  found  for  the 
more  distant  events.  Figure  1  shows  the  peak  frequency  versus  distance  for  all  three 
regional  phases.  This  trend  in  peak  frequency  versus  distance  is  consistent  enough  to  note 
even  though  events  of  the  same  magnitude  and  distance  can  have  significantly  different 
Pn  peak  frequencies.  For  example,  two  explosions  of  magnitude  2.6  at  a  distance  of  8 
degrees  can  have  peak  Pn  frequencies  of  2.0  or  10.0  Hz.  Whether  this  is  a  source  effect  or 
a  path  effect  is  yet  to  be  determined.  Regional  dependencies  certainly  exist.  Earthquakes 
on  the  west  coast  of  Norway  have  high  Pn  frequencies  (~10  Hz)  whereas  similar  sized 
events  at  the  same  distance  in  Sweden  are  of  lower  frequency  (~4  Hz). 

Q  estimates  for  Pn  are  generally  high,  with  no  obviously  azimuthal  dependence. 
Refinement  of  the  estimation  procedure  and  the  addition  of  more  events  should  resolve 
any  existing  trends  as  yet  unseen. 

The  Pn  waveforms  for  some  explosions  are  highly  monochromatic,  arriving  as  a 
wave  packet.  This  may  be  caused  by  delayed  firing  of  charges,  resulting  in  constructive 
interference  at  specific  frequencies.  Cepstral  analysis  of  these  waveforms  suggest  the  pres¬ 
ence  of  multiple  sources  separated  by  approximately  0.3  seconds.  The  presence  of  multi¬ 
ple  cepstral  peaks  shows  promise  as  a  means  of  identifying  mine  blasts,  however,  the 
interference  pattern  of  peaks  and  troughs  seen  in  the  spectra  of  these  events  can  severely 
bias  estimates  of  the  corner  frequency  and  Q. 

Pn  particle  motion  tends  to  be  very  rectilinear  for  a  few  seconds  immediately  follow¬ 
ing  the  initial  onset.  Pn  polarization  is  stronger  at  higher  frequencies  (around  10  Hz) 
whenever  there  is  a  good  signal-to-noise  ratio  in  this  band.  The  polarization  structure  of 
the  coda  waves  following  the  initial  Pn  arrival  deteriorates  in  time  so  that  the  late  P- 
coda  signals  just  prior  to  the  Sn  arrival  have  particle  motions  essentially  like  that  of 
ambient  noise.  This  P-coda  polarization  degrades  more  rapidly  at  higher  frequencies,  and 
the  duration  of  the  initial  Pn  strong  rectilinear  motion  decreases  with  increasing  source 
distance.  If  the  events  are  within  about  500  km  then  a  series  of  strongly-polarized  arrivals 
may  be  present  in  the  early  P- wave  coda  within  20  seconds  of  Pn,  corresponding  to  a 
series  of  P  waves  trapped  in  the  crust. 
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Polarization  attributes  of  the  Pn  motions  were  extracted  at  the  time  the  signal 
amplitude  and  rectilinearity  were  highest.  For  a  displacement  signal- to-noise  ratio  (S/N) 
greater  than  about  2,  PH  can  be  confidently  distinguished  from  noise  and  Lg  waves  on  the 
basis  of  rectilinearity.  However,  Sn  also  tends  to  be  somewhat  rectilinear  and  may  be  con¬ 
fused  with  Pn  on  this  basis,  depending  on  the  signal-to-noise  ratio.  The  best  way  to 
separate  Pn  from  noise,  Sm  and  Lf  appears  to  be  simply  the  vertical-to-horizontal  ratio 
(V/H)  at  the  highest  frequencies  for  which  the  signal-to-noise  ratio  of  Pn  is  greater  than 
about  1.5.  Figure  £  shows  a  plot  of  log(V/H)  versus  log(S/N)  for  about  40  regional 
events  recorded  at  NORESS.  The  signal-to-noise  values  in  Figure  £  are  in  amplitude 
units  and  the  vertical-to-horizontal  ratios  are  in  units  of  power.  Here  the  parameter  V /H 
is  computed  as  2  V2  /  (IT+N2),  where  V,  E  and  N  are  the  vertical,  east  and  north  ampli¬ 
tudes,  respectively.  Using  the  orientation  of  Pn  rectilinear  motion,  source  back  azimuths 
can  be  estimated  to  better  than  ±  20  *  in  about  90%  of  the  cases  when  the  S/N  ampli¬ 
tude  is  greater  than  about  1.5. 

S.1.4.  Sn 

Peak  frequencies  of  Sn  waves  were  in  all  cases  less  than  those  of  Pn  and  greater  than 
those  of  Lg.  These  ranged  from  1.5  Hz  to  9  Hz.  As  with  the  Pn  waves,  systematic  varia¬ 
tion  with  magnitude  did  not  show  up  in  our  dataset  although  there  is  a  weak  dependence 
on  distance  as  seen  in  Figure  1.  The  characteristics  of  the  Sn  waves  generally  followed 
those  of  their  associated  P n.  For  example,  when  the  Pn  peak  frequency  was  high,  the  Sn 
peak  frequency  was  similarly  high.  ThuB,  regional  dependencies  hold  for  both  phases.  It 
is  interesting  to  note  that  explosions  with  highly  monochromatic  Pn  waves  had  Sn  waves 
with  very  similar  spectra.  Considering  that  Pn  and  Sn  follow  nearly  identical  wave  paths 
in  the  Earth,  this  suggests  that  source  effects  may  play  a  prominent  role  in  the  character 
of  regional  phases.  Q  estimates  for  Sn  are  consistently  high  as  with  Pn  but  are  not  dis- 
cernibly  different  from  those  for  Pn  or  Lf  due  to  the  high  variance  associated  with  the 
estimates. 

Like  Pn,  the  Sn  particle  motions  are  strongly  polarized  over  a  short-time  interval  of 
only  a  few  seconds.  Sn  polarization  is  generally  not  as  strong  or  well-defined  as  Pn  polari¬ 
zation.  However,  the  analysis  carried  out  so  far  indicates  that  Sn  can  be  distinguished 
from  the  other  regional  phases  including  Lf  with  a  high  degree  of  confidence  using 
particle-motion  information.  The  Sn  phase  is  more  emergent  than  Pn,  and  its  signal-to- 
noise  level  is  often  lower  than  that  for  Pn>  since  it  is  the  late  P  coda  waves  which  control 
the  effective  background  noise  for  Sn.  Sn  polarization  is  generally  strongest  at  or  just 
prior  to  the  peak  three-component  envelope  of  this  phase.  Generally,  Sn  particle  motions 
are  somewhat  better  defined  at  the  higher  frequencies  than  at  the  lower  frequencies  (i.e., 
better  at  10  Hz  versus  4  Hz). 

The  dominant  particle-motion  characteristic  observed  for  Sn  is  the  large  horizontal 
motion  relative  to  vertical.  Figure  £  shows  the  vertical-to-horizontal  ratios  measured  for 
Sn  to  be  well  below  1.0  and  a  function  of  the  signal-to-noise  level.  The  scatter  is  quite 
large,  and  the  vertical-to-horizontal  ratio  in  itself  does  not  reliably  distinguish  S  from 
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Figure  2.  Vertical- to- horizontal  ratios  of  regional  phases  computed  from  three- 
component  seismograms.  Results  are  plotted  in  three  separate  frequency  bands  as  a  func¬ 
tion  of  the  signal-to-noise  ratio  in  each  band.  Pn  waves  tend  to  have  V/H  ratios  greater 
than  1.0,  Sn  waves  have  V/H  ratios  less  than  1.0,  and  waves  are  similar  to  the 
ambient  noise  with  a  V/H  near  1.0. 
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The  radial  motions  for  Sn  are  generally  about  equal  to  or  slightly  less  than  the  horizon¬ 
tals.  This  means  that  Sn  in  not  very  useful  for  estimating  source  azimuth.  The  near¬ 
equal  horizontals  and  smaller  verticals  gives  Sn  a  polarisation  structure  which  is  planar 
by  nature,  with  this  plane  of  motion  oriented  close  to  horizontal.  This  implies  that  the 
steeply-incident  Sn  waves  contain  both  and  Sy  components,  with  the  Sg  or 
transverse  component  somewhat  larger  in  many  cases. 

3.1.5.  Lg 

Peak  frequencies  of  L  waves  ranged  from  a  low  of  1.5  Hz  to  a  high  of  about  6  Hz, 
however,  the  cases  with  high  frequencies  were  very  limited.  Since  L  consists  predom¬ 
inantly  of  S-waves  traveling  and  reflected  within  the  crust,  they  are  highly  path  sensitive 
as  compared  with  Pn  and  Sn.  A  separation  in  peak  frequency  between  explosions  and 
microearthquakes  of  the  same  magnitude  did  not  show  up.  As  for  the  preceding  phases, 
Q  estimates  are  high,  and  peak  frequencies  show  a  dependence  on  distance  ( Figure  1). 

The  Lg  wavetrain  is  generally  longer  and  more  complex  than  the  other  phases,  with 
a  gradual  increase  in  the  envelope  before  the  peak  amplitude  is  reached.  particle- 
motions  are  quite  variable  with  time,  and  the  most  consistent  motions  are  found  around 
the  time  of  the  envelope  maximum.  Since  this  phase  is  longer  duration,  a  longer  window 
was  used  to  extract  Lg  particle-motion  information  than  was  used  for  Pn  or  Sn.  Figure  2 
shows  that  the  vertical  and  horizontal  amplitudes  of  L  tend  to  be  about  equal,  which  is 
also  true  of  ambient  noise.  However,  a  characteristic  feature  of  Lg  motions  is  the  domi¬ 
nance  of  vertical  and  transverse  components  relative  to  the  radial.  This  means  that,  like 
Sn,  the  Lg  polarization  structure  is  relatively  planar.  Unlike  Sn,  the  orientation  of  this 
plane  of  motion  is  dominantly  vertical  and  transverse  rather  than  horizontal.  The 
analysis  carried  out  so  far  indicates  Lg  can  be  separated  from  Sn  on  the  basis  of  the  orien¬ 
tation  of  this  plane  of  polarization,  or,  equivalently,  on  the  basis  of  the  radial-to-vertical 
amplitude  ratios.  Such  observed  Lg  polarization  structure  is  consistent  with  the  theory 
that  the  Lg  wavetrain  in  the  short-period  band  is  comprised  of  trapped  S  waves  propagat¬ 
ing  relatively  horizontally  in  the  crust  with  both  Sy  and  SH  present. 

The  dominance  of  transverse  motion  for  Lg  implies  that  the  source  azimuth  can  be 
estimated  using  the  orientation  of  minimum  horizontal  motion.  The  minimum  axis  of  the 
polarization  ellipse  computed  from  only  the  horizontal  components  has  been  previously 
used  to  estimate  the  source  back  azimuth  (to  within  ±180*),  (Smart  and  Sproules,  1981). 
Preliminary  analysis  of  regional  events  at  NORESS  indicates  that  the  source  azimuth  can 
be  determined  this  way  to  within  about  20  *  in  more  than  75%  of  the  cases  when  the  S/N 
of  Lg  is  greater  than  1.5.  Thus,  the  source  azimuth  estimates  from  Lg  are  not  as  accurate 
as  those  obtained  from  the  P_  motions. 


A.  Jurkevics,  P.  Dysart  and  J.  J.  Pulli 
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4.  INTERNATIONAL  DATA  EXCHANGE 


4.1.  GSE  WAVEFORM  EXCHANGE  EXPERIMENTS 


4.1.1.  Introduction 

Recent  activities  of  the  GSE  in  Geneva  have  accelerated  the  development  of  capabil¬ 
ities  to  exchange  seismic  waveforms.  A  data  format  tentatively  approved  by  the  GSE  in 
August  1986  was  reviewed  and  revised  during  the  latter  part  of  1986  and  the  final  version 
was  distributed  to  GSE  members  in  late  December  1986. 

This  report  documents  waveform  exchange  experiments  and  related  activities  con¬ 
ducted  by  the  Center  during  the  first  quarter  of  1987  in  order  to  gain  experience  with  the 
proposed  procedures  and  formats,  and  provides  observations  of  Center  personnel  on  the 
various  exchange  procedures. 

4.1.2.  Waveform  Data  Exchange  Methods 

The  Center  has  conducted  experiments  employing  all  of  the  digital  waveform 
exchange  methods  proposed  by  the  GSE.  These  have  included  transfers  by  magnetic  tape, 
electronic  mail  and  direct  computer-to-computer  file  transfer.  The  waveform  data 
transfers  are  summarized  in  Table  I.  Each  exchange  method  has  its  advantages  and 
disadvantages,  and  these  are  discussed  briefly  here.  It  should  be  noted  that  the  way  in 
which  a  transfer  is  carried  out  can  have  profound  implications  for  cost,  reliability,  data 
quality,  and  handling  of  the  received  data  files. 

The  WMO/GTS  has  been  used  for  years  to  transfer  amounts  of  data  that  are  small 
in  comparison  to  useful  seismic  waveforms,  and  it  was  used  exclusively  for  the  exchange  of 
parameter  data  during  the  GSETT  experiment  in  1984.  It  is  straightforward  to  broadcast 
copies  of  data  to  many  nations  over  the  WMO/GTS  network.  There  are  doubts  that  the 
WMO/GTS  will  be  able  to  handle  the  volume  of  data  that  would  be  expected  from  a 
large-scale  exchange  of  waveforms.  Furthermore,  since  the  WMO  system  does  not  have 
error  correction  capability,  and  waveform  data  is  useless  unless  almost  error  free,  data 
integrity  could  be  a  serious  problem. 

Use  of  electronic  mail  (El-mail)  networks  has  been  limited  by  computer  compatabil- 
ity  and  the  coverage  of  the  individual  networks.  El-mail  can  be  transferred  over  standard 
telephone  lines  or  packet-switched  networks,  and  it  works  as  long  as  there  is  a  path  and 
address  available  to  the  sender.  Use  of  El-mail  systems  can  reduce  costs  by  permitting 
scheduling  of  transfers  to  take  advantage  of  periods  of  low  transmission  rates  and  by 
spooling  transfers.  Electronic  mail  systems  can  also  incorporate  data  compression  and 
quality  assurance  algorithms. 


INTERNATIONAL  WAVEFORM  EXCHANGE  EXPERIMENTS 
WITH  UNITED  STATES  PARTICIPATION 
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Direct  computer-to-computer  transfers  (excluding  E-mail)  are  usually  initiated  by  an 
individual  while  logged  onto  a  remote  computer.  This  method  of  transfer  can  be  wasteful 
for  a  number  of  reasons.  First,  every  user  must  have  login  ca^  -iOility  and,  depending  on 
transfer  direction,  write  permission  on  the  remote  machine.  Second,  transfers  of  this  type, 
once  initiated,  tie  up  terminals,  ports  and  people.  Computer-to-computer  connections 
designed  for  human  interaction  frequently  Me  set  up  to  pass  each  character  as  it  is  typed. 
This  can  cause  packet-switched  network  charges  to  be  orders  of  magnitude  higher  than 
they  might  otherwise  be  if  the  connection  was  optimized  for  transfer  without  interaction. 
Third,  unless  files  are  kept  small  or  another  layer  of  protocols  are  used  (e.g.,  KERMIT), 
retransmission  costs  can  be  significant.  Thus,  this  method  of  transfer  can  increase  the 
data  auditing  task  associated  with  exchanges.  Since  this  method  allows  a  user  to  initiate 
a  transfer  to  his  own  machine  while  logged  in  on  a  remote  machine,  it  introduces  an  addi¬ 
tional  class  of  transfer  which  contributes  to  an  existing  problem  concerning  all  transfers. 
That  is,  having  the  final  destination  and  name  of  the  data  be  independent  of  the  transmis¬ 
sion  method.  A  solution  to  this  problem  is  outlined  in  a  subsequent  section.  Finally,  mag¬ 
netic  tape,  although  a  delayed  method  for  transferring  data,  is  widely  used  and  highly 
reliable  and  does  not  require  special  communications  facilities  or  on-line  storage  which 
may  be  a  problem  for  some  centers.  The  cost  of  tape  transfer  is  competitive  with 
machine-to-machine  transfers. 


4.1.3.  Data  Requests  and  Exchanges 


Although  it  is  expected  that  the  exchange  of  waveforms  will  be  initiated  by  a  formal 
data  request,  and  the  GSE  has  specified  a  format  for  data  requests  as  well  as  for  the  data 
itself,  most  of  the  trial  exchanges  were  not  initiated  with  a  formal  request.  The  Center 
received  only  one  request  in  GSE  format  and  that  came  from  Japan.  The  request  was 
incomplete  and  a  clarification  was  requested.  The  clarification  was  also  incomplete  since 
the  time  given  was  arrival  time  rather  than  start  time  and  no  window  lengths,  or  phase 
names  provided.  But  some  assumptions  were  made  and  the  requested  data  were  E-mailed 
to  Japan’s  account  on  seismo  within  hours.  Similarly,  data  transfers  to  the  Center 
resulted  from  a  non-standard  blanket  request  for  waveform  data  from  three  widely 
recorded  earthquakes.  The  request  was  forwarded  along  with  tapes  containing  waveform 
data  in  GSE  format,  which  were  distributed  to  all  GSE  members  in  January  1987  for  for¬ 
mat  evaluation.  It  was  this  blanket  request  which  produced  the  exchanges  summarized 
here. 


Table  I  summarizes  waveform  exchanges  which  occurred  during  the  first  quarter  of 
1987.  Information  on  the  quantity  of  data  transmitted,  the  transmission  rate  where  it 
could  be  determined,  method  of  transmission  and  event(s)  for  which  data  were  transmit¬ 
ted  are  included  in  the  summary.  The  events  are  designated  by  the  year  and  day  of  year. 
Table  II  gives  detailed  information  on  the  events  selected  for  waveform  data  exchange. 
The  Center  received  five  GSE  format  exchange  tapes.  A  non-standard  format  tape  was 
also  received  from  Czechoslovakia.  The  GSE  tapes  and  their  contents  are  tabulated  in 
Table  III.  The  entries  are  images  of  the  GSE  standard  waveform  identification  header 
which  begins  with  the  character  string  WIDn  and  gives  the  event  date,  start  time, 
number  of  samples,  station  name,  channel  identifiers,  sample  rate,  instrument  type  and 
data  format  type.  Waveform  files  that  have  been  delivered,  retrieved  or  E-mailed  to 
seismo  are  tabulated  in  Table  IV. 
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Table  HI 


CONTENTS  OF  GSE  TAPES  RECEIVED: 
AUSTRIA: 


WID1 

1984338 

4 

20 

12 

800 

2400 

KBA 

SP2 

SZ 

20.0000000 

S-13 

INT8 

CANADA: 

WID1 

1985009 

10 

34 

15. 

449 

401 

YKA 

SZ 

20.0000000 

INT10 

WID1 

1984345 

10 

31 

34. 

966 

3001 

GAC 

SZ 

30.0000000 

INT10 

WIDl 

1984345 

10 

31 

34. 

966 

3001 

GACSP 

SN 

30.0000000 

INT10 

WID1 

1984345 

10 

31 

34. 

966 

3001 

GACSP 

SE 

30.0000000 

INT10 

WIDl 

1984345 

10 

19 

36.000 

1801 

GAC 

LZ 

0.5000000 

INTlO 

WIDl 

1984345 

10 

19 

36.000 

1801 

GAC 

LN 

0.5000000 

INT10 

WIDl 

1984345 

10 

19 

36.000 

1801 

GAC 

LE 

0.5000000 

INTlO 

ITALY 

WIDl 

1984342 

10: 

28 : 

51. 

.880 

7203 

MNS 

SPZ 

SZ 

50.0000000 

S-13 

INT5 

JAPAN 

WIDl 

1984338 

13 

10 

16 

960 

9000 

MAT 

SPZ 

SZ 

20.0000000 

JAPAN 

INT 

WIDl 

1984338 

13 

10 

16 

960 

9000 

MAT 

SPN 

SN 

20.0000000 

JAPAN 

INT 

WIDl 

1984338 

13 

10 

16 

960 

9000 

MAT 

SPE 

SE 

20.0000000 

JAPAN 

INT 

WIDl 

1984338 

13 

4 

58 

960 

8850 

MAT 

LPZ 

LZ 

1.0000000 

JAPAN 

INT 

WIDl 

1984338 

13 

4 

58 

960 

8850 

MAT 

LPN 

LN 

1.0000000 

JAPAN 

INT 

WIDl 

1984338 

13 

4 

58 

960 

8850 

MAT 

LPE 

LE 

1.0000000 

JAPAN 

INT 

WIDl 

1984339 

16 

46 

1 

960 

9000 

MAT 

SPZ 

SZ 

20.0000000 

JAPAN 

INT 

WIDl 

1984339 

16 

46 

1 

960 

9000 

MAT 

SPN 

SN 

20.0000000 

JAPAN 

INT 

WIDl 

1984339 

16 

46 

1 

960 

9000 

MAT 

SPE 

SE 

20.0000000 

JAPAN 

INT 

WIDl 

1984339 

16 

40 

40 

960 

8850 

MAT 

LPZ 

LZ 

1.0000000 

JAPAN 

INT 

WIDl 

1984339 

16 

40 

40 

960 

8850 

MAT 

LPN 

LN 

1.0000000 

JAPAN 

INT 

WIDl 

1984339 

16 

40 

40 

960 

8850 

MAT 

LPE 

LE 

1.0000000 

JAPAN 

INT 

WIDl 

1984345 

19 

40 

47 

960 

9000 

MAT 

SPZ 

SZ 

20.0000000 

JAPAN 

INT 

WIDl 

1984345 

19 

40 

47 

960 

9000 

MAT 

SPN 

SN 

20.0000000 

JAPAN 

INT 

WIDl 

1984345 

19 

40 

47 

960 

9000 

MAT 

SPE 

SE 

20.0000000 

JAPAN 

INT 

SWEDEN: 

WIDl 

1984338 

04 

19 

20 

000 

4070 

HFS 

HFS 

SZ 

60.0000000 

S13 

INT6 

WIDl 

1984338 

04 

19 

20 

000 

4070 

HFS 

APP 

SZ 

60.0000000 

S13 

INT6 

WIDl 

1984338 

0 

19 

20 

000 

4070 

HFS 

SLL 

SZ 

60.0000000 

S13 

INT6 

WIDl 

1984338 

04 

19 

20 

000 

4070 

SKI 

KTB 

SZ 

60.0000000 

S13 

INT6 

WIDl 

1984338 

04 

19 

20 

000 

4070 

SKI 

HRN 

SZ 

60.0000000 

S13 

INT6 

WIDl 

1984339 

07 

55 

50 

000 

4070 

HFS 

HFS 

SZ 

20.0000000 

S13 

INT6 

WIDl 

1984339 

07 

55 

50 

000 

4070 

HFS 

SLL 

SZ 

20.0000000 

S13 

INT6 

WIDl 

1984339 

07 

55 

50 

000 

4070 

HFS 

APP 

SZ 

20.0000000 

S13 

INT6 

WIDl 

1984339 

07 

55 

50 

000 

4070 

HFS 

TBY 

SZ 

20.0000000 

S13 

INT6 

WIDl 

1984338 

4 

20 

27 

927 

3599 

RSNY 

SPZ 

SZ 

39.9998055 

S-750 

INTV 

WIDl 

1984338 

4 

30 

13 

192 

3599 

RSNY 

LPZ 

LZ 

0.9999952 

KS-360 

INTV 

WIDl 

1984338 

4 

30 

13 

192 

3599 

RSNY 

LPN 

LN 

0.9999952 

KS-360 

INTV 

WIDl 

1984338 

4 

30 

13 

192 

3599 

RSNY 

LPE 

LE 

0.9999952 

KS-360 

INTV 

WIDl 

1984147 

12 

20 

37 

960 

3900 

MAT 

SPZ 

SZ 

20.0000000 

JAPAN 

INTV 

WIDl 

1984147 

12 

15 

31 

960 

8850 

MAT 

LPZ 

LZ 

1.0000000 

JAPAN 

INTV 

WIDl 

1984147 

12 

15 

31 

960 

8850 

MAT 

LPN 

LN 

1.0000000 

JAPAN 

INTV 

WIDl 

1984147 

12 

15 

31 

960 

8850 

MAT 

LPE 

LE 

1.0000000 

JAPAN 

INTV 

WIDl 

1984345 

10 

31 

34 

966 

3000 

GAC 

SPZ 

SZ 

30.0000000 

KS-36K 

INTlO 

WIDl 

1984345 

10 

19 

36 

000 

1728 

GAC 

LPZ 

LZ 

0.5000000 

GS-36K 

INTlO 

WIDl 

1984345 

10 

19 

36 

000 

1800 

GAC 

LPN 

LN 

0.5000000 

KS-36K 

INTlO 

WIDl 

1984345 

10 

19 

36 

000 

1800 

GAC 

LPE 

LE 

0.5000000 

KS-36K 

INTlO 

WIDl 

1985009 

10 

34 

15 

449 

400 

YKA 

SPZ 

SZ 

20.0000000 

MARK I I 

INTlO 

WIDl 

1985355 

10 

3 

27 

104 

5865 

CLL 

SPZ 

SZ 

19.5300007 

CIPE 

INTV 

WIDl 

1986249 

23 

15 

12 

0 

3964 

WMQ 

SPZ 

SZ 

40.0000000 

STSl 

INTV 
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Table  IV 


FILES  ACQUIRED  DURING  FIRST  QUARTER  1987: 

file  name  how  acquired 


/a/gse/denmark/wexltestl 

del 

i  vered 

/a/gse/denmark/wexltest2 

del 

i  vered 

/a/gse/japan/84123s.mat 

del 

i vered 

/a/gse/japan/84147L.MAT 

del 

i vered 

/a/gse/japan/84147s.mat 

del 

i vered 

/ a/gse/ j apan/84152s . mat 

del 

i vered 

/a/gse/japan/84196s.mat 

del 

i vered 

/a/gse/japan/843381 .mat 

del 

i vered 

/a/gse/japan/84338s.mat 

del 

i vered 

/a/gse/japan/843391 .mat 

del 

i vered 

/a/gse/ j  apa n/84339s . ma  t 

del 

i vered 

/a/gse/japan/84345s.mat 

del 

i vered 

/a/gse/norway / 1  eve  1 2 . samp  1 e 

del 

i vered 

/a/gse/ sweden/ev841203 . dat 

del 

i vered 

/a/gse/sweden/ev841204.swe 

del 

i vered 

/a/gse/uk/testl 

del 

i vered 

/a/gse/uk/ukdec03.gse 

del 

i vered 

/a/gse/uk/ukdec04.gse 

del 

i vered 

/a/gse/uk/ukdeclO . gse 

del 

i vered 

/a/gse/usa/f rgerman.data 

Ema 

i  led 

/a/gse/usa/gac I e. tk 
/a/gse/usa/gac I n . tk 
/a/gse/usa/gacsz . tk 
/a/gse/usa/ykasz . tk 


retrieved 
retr i eved 
retrieved 
retrieved 


This  table  shows  one  potential  problem  with  the  uncontrolled  exchanges  which  have 
been  conducted  thus  far.  Those  files  delivered  by  nations  exist  in  their  home  directory 
and  are  given  names  by  the  delivering  nation.  Files  E-mailed  or  retrieved  are  named  by 
the  Center  staff  member  collecting  them,  but  cannot  be  placed  in  the  directory  of  the 
country  of  origin  (no  write  permission),  so  they  are  placed  in  the  home  directory  of  the 
U.S.A.  account.  What  should  be  developed  is  a  system  for  naming  the  waveform  files,  as 
well  as  the  designation  of  a  waveform  storage  area  with  a  directory  for  each  nation  where 
files  can  be  deposited  by  Center  staff.  After  a  nation  delivers  a  file,  the  Center  should  be 
notified  by  E-mail  so  that  the  file  can  be  moved  from  the  sending  nation’s  home  directory 
into  the  waveform  storage  area.  For  the  situation  where  the  file  arrives  via  E-mail  or  is 
retrieved  by  the  Center,  the  Center  staff  names  the  file  and  places  it  in  the  sending 
nation’s  directory  in  the  waveform  area. 

4.1.4.  Problems/Solutions 

Experience  has  shown  that  the  default  PAD  (packet  assembler  disassembler)  settings 
frequently  are  not  optimum.  Nations  wishing  to  transmit  data  should  check  for  the 
correct  PAD  settings  prior  to  each  data  transmission  to  insure  optimum  transfer  rates. 
Incorrect  PAD  settings  may  increase  data  transmission  cost  by  a  factor  of  50  or  more 
over  costs  when  using  correct  settings.  Moreover,  the  settings  for  one  computer  may  not 
be  optimum  for  other  computers,  suggesting  that  records  should  be  kept  of  optimum  PAD 
settings  for  each  transmission  path. 

To  maintain  an  adequate  record  of  data  exchanges,  it  is  important  that  unique  iden¬ 
tification  codes  be  incorporated  both  in  requests  for  waveforms  and  in  the  waveform  data 
files  themselves.  The  unique  request  identifier  should  be  assigned  by  the  nation  originat¬ 
ing  a  request,  and  the  data  identifier  should  be  assigned  by  the  nation  providing 
waveform  data.  Each  nation  should  be  responsible  to  insure  that  the  numbers  it  sends 
out  are  unique.  The  GSE  should  establish  a  coding  system  which  can  be  used  by  the 
NDCs  and  IDCs  to  uniquely  identify  waveform  requests  and  waveform  transmissions.  It 
is  suggested  here  that  the  data  id  part  of  the  WEX  header  become  the  file  name  of  that 
data  with  wex  prepended.  Furthermore,  it  is  suggested  that  the  id  incorporate  the  three- 
character  ISO  codes  given  in  Table  V.  Thus,  for  a  data  package  (wex,  cal,  wid  and  dat 
sections)  provided  by  the  U.S.A.  with  the  data  id  USA001290  would  have  the  disc- 
resident  file  name  WEXUSA001290.  This  naming  convention  could  also  extend  to  physi¬ 
cal  tape  labels. 

It  is  further  suggested  that  the  GSE  establish  a  procedure  for  submitting  requests 
and  a  place  for  depositing  data  on  each  machine  assigned  the  GSE  duties.  When  data  is 
deposited,  a  notice  should  be  sent  to  the  receiving  nation  via  E-mail.  It  is  suggested  that 
nations’  accounts  and  directory  names  also  be  based  on  the  ISO  three-letter  codes  given  in 
Table  V  rather  than  on  the  full  name  of  the  nation  or  some  ad  hoc  abbreviation. 
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Table  V 


NATION 

ISO 

2char 

CODE 

3char 

AUSTRALIA 

AU 

AUS 

AUSTRIA 

AT 

AUT 

BELGIUM 

BE 

BEL 

BULGARIA 

BG 

BGR 

CANADA 

CA 

CAN 

CZECHOSLOVAKIA 

CS 

CSK 

DENMARK 

DK 

DMK 

EGYPT 

EG 

EGY 

FINLAND 

FI 

FIN 

GERMAN  DEMOCRATIC  REPUBLIC 

DD 

DDR 

FEDERAL  REPUBLIC  OF  GERMANY 

DE 

DEU 

HUNGARY 

HU 

HUN 

INDIA 

IN 

IND 

INDONESIA 

ID 

IDN 

IRAN 

IR 

IRN 

ITALY 

IT 

ITA 

JAPAN 

JP 

JPN 

NETHERLANDS 

NL 

NLD 

NEW  ZEALAND 

NZ 

NZL 

NORWAY 

NO 

NOR 

PERU 

PE 

PER 

ROMANIA 

RO 

ROM 

SWEDEN 

SE 

SWE 

UNITED  KINGDOM 

GB 

GBR 

UNITED  STATES 

US 

USA 

USSR 

SU 

SUN 

ZAMBIA 

ZM 

ZMB 
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4.1.5.  Other  Experiment! 

One  ancillary  experiment  was  conducted  in  order  to  determine  how  many  station- 
events  can  be  written  to  a  single  tape.  A  standard  2400-foot  reel  was  used  and  the  tape 
was  written  at  1600  bpi.  A  station-event  is  defined  here  as  a  set  of  four  waveforms  (sz 
and  3  lp)  in  standard  GSE  length,  each  preceded  by  its  own  calibration  section.  One 
station-event  is  written  per  tape  file.  Data  from  a  single  RSTN  station  was  repeatably 
written  with  a  complete  set  of  four  waveforms  (sz,  lz,  In,  le)  and  a  calibration  section 
included  in  each  file.  Waveforms  were  in  the  standard  GSE  length  (90  sec  sz,  3600  sec 
lp).  288  station-events  were  written  to  the  tape.  Given  a  50-station  network,  it  is  then 
reasonable  to  assume  that  one  could  write  data  for  five  events  onto  a  similar  tape  (250 
station-events).  This  assumes  that  the  calibration  sections  from  the  50  stations  would  be 
of  a  similar  size  to  the  one  for  a  RSTN  station.  The  implication  here  is  that  if  all  50  sta¬ 
tions  exchanged  data  for  40  events  a  day,  eight  tapes  a  day  would  be  required  to  store 
such  data.  This  probably  represents  an  upper  limit  since  an  average  of  20  teleseisms  were 
recorded  each  day  during  the  GSETT  and  not  all  of  these  events  were  recorded  at  all  sta¬ 
tions.  At  the  other  extreme,  if  five  events  per  week  were  selected  for  exchange,  one  tape 
per  week  would  be  required. 

4.1.6.  New  Software 

Major  modifications  to  the  GSE  Level  II  file  reading  program  ‘r2’  were  made  during 
the  last  quarter.  While  no  major  changes  to  the  format  precipitated  these  changes,  they 
are  worth  mentioning.  r2  now  has  the  ability  to  check  data  files  for  compliance  with  GSE 
rules.  While  the  check  is  far  from  exhaustive,  it  will  prevent  random  files  from  being  pro¬ 
cessed  and  producing  erroneous  databases.  Data  files  are  tested  to  insure  that  the  key¬ 
words  ‘WEXl’,  ‘CALI’,  ‘WID1’,  and  ‘DATl’  are  in  place.  The  number  of  data  values  in 
the  DAT  section  are  checked  to  insure  they  equal  the  number  of  samples  given  in  the 
WID  section.  While  it  is  not  a  hard  error  for  a  file  to  lack  a  WEX  or  CAL  section,  files 
without  WID  sections  are  considered  non-standard.  A  manual  page  has  been  prepared  for 
r2,  but  the  software  has  not  been  officially  installed  since  it  is  expected  that  further 
modifications  will  be  needed  as  we  gain  experience  with  the  format. 


Michael  Tiberio 
George  Bullin 
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5.  OTHER  RESEARCH,  STATUS  OF  DATABASES,  SOFTWARE 
DEVELOPMENT 


5.1.  SEN SITIVITIY  OF  mh  TO  EARTH  ATTENUATION 


5.1.1.  Outline 

Work  is  currently  in  progress  on  a  study  to  investigate  the  sensitivity  of  mft  to  earth 
attenuation  using  teleseismic  recordings  of  presumed  nuclear  explosions.  This  is  a  follow¬ 
up  to  a  previous  study  (Jurkevics  and  Romney,  1986)  in  which  seismograms  were  simu¬ 
lated  by  combining  a  source  model,  near-source  reflection  coefficients,  earth  attenuation 
and  an  instrument  response.  The  current  study  involves  real  seismograms  recorded  glo¬ 
bally  from  presumed  nuclear  explosions  at  various  sites  in  the  U.S.  and  the  U.S.S.R.  The 
basic  idea  is  to  take  the  recorded  seismograms,  apply  additional  earth  attenuation  by  a 
convolution  operation  and  to  then  measure  the  change  in  amplitude,  period  and  associ¬ 
ated  my  This  approach  of  adding  additional  attenuation  to  recorded  signals  is  justified 
by  the  linearity  of  the  convolution  operator  with  which  the  earth  attenuation  is  applied. 
Several  attenuation  models  are  used  including  constant  Q  and  frequency-dependent  Q, 
the  latter  based  on  a  linear  attenuation  model  proposed  by  Minster  (1978). 

Short-period  recordings  of  a  total  of  42  events  from  12  different  geographic  areas  are 
being  analyzed.  The  events  include  U.S.  nuclear  tests  and  presumed  U.S.S.R.  weapons 
tests  and  peaceful  nuclear  explosions  (PNEs).  The  source  areas  include  the  Nevada  Test 
Sites  (NTS),  Amchitka  Island,  Degelen  Mt.  and  Shagan  River  sites  in  East  Kazakh, 
Novaya  Zemlya,  Azgir,  Lake  Baykal  region,  Kushata,  and  a  few  others.  The  recordings 
were  obtained  from  several  seismic  networks  including  WWSSN,  DWWSSN,  SRO,  ASRO, 
and  RSTN.  The  stations  were  selected  to  provide  representative  global  coverage,  con¬ 
strained  by  the  available  data  in  the  database  at  the  Center  for  Seismic  Studies.  Usually 
about  15  recordings  were  analyzed  for  each  event.  This  note  presents  a  few  initial  results 
and  illustrs  *  ss  the  nature  of  the  study;  a  more  detailed  account  will  be  presented  in  a  sub¬ 
sequent  report. 

5.1.2.  Preliminary  Results 

Figure  1  illustrates  the  change  in  for  different  amounts  of  attenuation  and  for 
three  different  source  sizes,  as  represented  by  the  three  different  corner  frequencies.  This 
result  is  for  simulated  seismograms  and  is  similar  to  the  results  presented  in  the  previ¬ 
ously  published  part  of  the  study  (Jurkevics  and  Romney,  1986).  The  source  model  is 
from  von  Seggern  and  Blandford  (1972),  with  an  overshoot  parameter  of  0.5.  A  surface 
reflection  coefficient  of  -0.5  at  a  delay  of  0.25  seconds,  a  constant  Q  model,  and  the 
WWSSN  short-period  instrument  response  were  used  in  this  particular  example.  The  mb 
values  in  Figure  1  have  been  normalized  to  0.0  at  t*  —  0.0  in  order  to  emphasize  the  sen¬ 
sitivity,  or  slope,  of  mt  relative  to  t*.  Figure  1  shows  that,  as  expected,  the  values 
from  sources  with  lower  corner  frequency  are  less  sensitive  to  added  absorption  (increased 
t*)  than  are  values  from  high  frequency  sources.  Further,  for  a  given  source,  the 


change  in 


WWSSN 


constant  Q 


Figure  1.  Sensitivity  of  to  attenuation  as  measured  from  synthetic  seismograms.  A 
von  Seggern-Blandford  source  pulse  using  three  different  corner  frequencies  was  tested. 
The  instrument  response  is  WWSSN  short-period  and  the  Q  model  is  independent  of  fre¬ 
quency.  The  m4  values  have  been  normalized  at  t  *  =  0  in  order  to  compare  the  slopes  of 
the  curves.  The  slope  is  a  function  of  source  corner  frequency  and  the  amount  of  attenua¬ 
tion  present. 


change  in  magnitude  (  Am(  )  from  additional  absorption  (  At*  )  is  a  function  of  t* 
itself.  Thus,  the  change  in  magnitude  resulting  from  transmission  of  P- waves  through 
some  particular  absorptive  layer  will  depend  on  both  the  spectrum  of  the  source  and  on 
the  absorption  along  the  remainder  of  the  path.  Changes  in  the  predominant  period  of  P 
also  result  from  changes  in  t  *,  as  noted  in  our  earlier  paper. 

We  apply  these  ideas  to  real  seismograms  from  explosions  as  a  partial  test  of 
assumptions  regarding  characteristics  of  sources  at  several  test  sites,  as  well  as  current 
concepts  of  absorption  under  those  test  sites.  For  instance,  if  we  assume  the  NTS  and 
Kazakh  sources  are  identical  (at  a  given  yield)  but  that  the  absorption  along  the  paths  to 
a  seismic  station  differs,  then  applying  an  additional  absorption  to  the  signals  should 
result  in  changes  in  that  differ  between  NTS  and  Kazakh.  Further,  it  should  be  possi¬ 
ble  to  transform  P  waveforms  from  one  site  to  model  the  main  characteristics  of  those 
from  the  other  by  application  of  more  attenuation. 

The  procedure  is  illustrated  in  Figure  2  which  shows  six  sets  of  waveforms  from  six 
different  events  recorded  at  one  station;  WWSSN  station  KTG  in  eastern  Greenland. 
Each  set  in  Figure  2  contains  three  signals;  the  top  trace  is  the  short-period  raw  seismo¬ 
gram  of  the  early  P- wave.  The  two  other  traces  in  each  set  are  the  same  signals  with  dif¬ 
ferent  amount  of  additional  attenuation  added.  In  this  case  the  additional  attenuation 
consists  of  a  constant  Q  model  with  a  t*  of  0.2  and  0.6  seconds.  The  three  sets  of 
seismograms  along  the  top  in  Figure  2  are  from  the  Nevada  Test  Site  with  reported  mag¬ 
nitudes  of  6.0,  5.5  and  5.6.  The  three  sets  below  are  presumed  tests  in  East  Kazakh  with 
magnitudes  of  6.0,  5.9  and  5.6.  The  peak  amplitude  in  nanometers  is  displayed  adjacent 
to  each  waveform.  The  main  effects  of  adding  the  additional  attenuation  are  the  decrease 
in  signal  amplitude  and  increase  of  low  frequency  energy  relative  to  high  frequency.  In 
the  analysis  the  values  were  measured  interactively  on  a  SUN  workstation.  The  so- 
called  "c"  measure  was  used,  which  is  the  amplitude  from  the  first  trough  to  the  subse¬ 
quent  peak.  While  data  from  a  single  station  are  insufficient  for  drawing  firm  conclu¬ 
sions,  these  figures  illustrate  the  previously  reported  result  that  NTS  explosions  (possibly 
aside  from  Piledriver)  differ  from  Kazakh  explosions  in  ways  that  are  not  explained  by 
greater  absorption  under  NTS. 

Figure  S  shows  preliminary  results  of  analyzing  actual  seismograms  for  eight 
presumed  PNEs  located  in  the  Azgir  region  of  U.S.S.R.  The  parameter  "slope"  in  the 
!  equation  A  mfc  =  slope  At*  is  plotted;  this  was  determined  by  adding  attenuation  of 

!  t*  =  0.2seconds  to  the  seismograms.  Both  constant  and  frequency-dependent  models  of 

|  Q  were  used.  The  recordings  were  from  RSTN  and  SRO  stations  with  an  instrument 

1  correction  applied  to  convert  to  WWSSN  short-period  response.  The  station  codes  are 

shown  on  the  graphs  and  the  magnitudes  of  the  events  are  shown  below.  The  values  plot¬ 
ted  for  each  event  in  Figure  2  show  a  considerable  scatter  between  stations.  However,  the 
trend  is  to  decrease  the  sensitivity  of  mt  to  added  t*  for  increasing  source  size;  this  is 
i  qualitatively  consistent  with  the  expected  decrease  in  corner  frequency  as  the  yield  (and 

|  magnitude)  increases. 


Andy  Jurkevics  and  Carl  Romney 


WWSSN  Station  KTG  (Eastern  Greenland) 


t*  added 
.0 

.2 

.6 

t*  added 
.0 

2 

.6 


75.134  NTS  6.0 


69.33 

~  A 

4  Shogon  6.0 

79.8 

I 

|/A/  Vv  W  v 

|  340 

1  A  A  J\m  Aa.  A 

“"I 

A 

L 

/ur 

A^l 

r 

1  \ /  V  V 

0  5  10 


70.287  NTS  5.5 


71.115  Degelen  5.9 
137. 


75.117  Shogon  5.6 


0  5 

time  (seconds) 


10 


Nevada 


East  Kazakh 


Figure  2.  Attenuation  is  added  to  short-period  P-waves  from  six  events  recorded  at 
WWSSN  station  KTG.  The  peak  amplitudes  in  nanometers  are  shown  adjacent  to  each 
trace.  The  year.day-of-year,  location,  and  m.  magnitude  are  given  for  each  event. 
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Figure  S.  Observed  values  for  "elope"  in  the  relation  Am^  =  elope ‘At*  are  shown  for 
eight  presumed  PNEs  in  the  Azgir  region  of  U.S.S.R.  The  recordings  are  from  the  RSTN 
and  SRO  networks.  Instrument  corrections  were  made  to  correct  to  WWSSN  short- 
period  response  first.  Both  constant  and  frequency-dependent  Q  models  were  used.  The 
year.day-of-year  and  mt  magnitudes  are  shown  for  each  event.  The  sensitivity  is 
larger  for  the  constant  Q  model  and  some  dependence  with  magnitude  is  also  apparent. 
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Dr.  Robert  Blandford  provided  valuable  comments  during  the  course  of  this  study. 
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5.2.  BODY  WAVE  MAGNITUDES  OF  EASTERN  KAZAKH 
EXPLOSIONS  CALCULATED  WITH  NORESS  DATA 


5.2.1.  Introduction 

Recently  the  Soviet  Union  ended  its  self-imposed  moratorium  on  nuclear  weapons 
testing.  Subsequently,  three  tests  were  conducted  at  its  Eastern  Kazakh  test  site,  approx¬ 
imately  38  degrees  from  the  NORESS  array.  The  Center  for  Seismic  Studies  receives  data 
from  the  array  in  real  time  via  satellite.  The  NORESS  data  were  analyzed  after  each 
test,  and  one  of  the  first  source  parameters  determined  was  the  body  wave  magnitude. 
The  mb s  calculated  from  the  NORESS  data  were  higher  than  those  calculated  from 
RSTN  stations  or  reported  by  the  U.S.  Geological  Survey’s  National  Earthquake  Informa¬ 
tion  Service  (NEIS).  The  path  between  Eastern  Kazakh  and  NORESS  is  an  efficient  pro¬ 
pagator  of  P- wave  energy  (Ringdal  and  Husebye,  1982).  However,  the  mb  residuals  for 
the  three  tests,  +0.38,  +0.74,  and  +1.1  mb  units,  were  considered  large  enough  and  dif¬ 
ferent  enough  to  warrant  further  review.  In  this  study,  we  calculated  body  wave  magni¬ 
tudes  for  ten  Eastern  Kazakh  explosions. 

5.2.2.  Data  and  Analysis 

The  ten  explosions  studied  are  listed  in  Table  I.  Epicentral  information  and  mb s  are 
taken  from  NEIS’s  Monthly  Summaries,  with  the  exception  of  the  1987  events  whose  epi¬ 
central  parameters  were  taken  from  NEIS’s  Quick  Epicenter  Determination  list.  Beamed 
P-waveforms  for  each  event  are  shown  in  Figure  1 ,  arranged  in  increasing  order  of  magni¬ 
tude.  The  waveforms  show  some  variation  with  magnitude,  which  may  reflect  the 
increasing  depth  of  burial.  Peak-to-peak  ground  motions  were  measured  from  the  second 
peak  and  second  trough  of  each  waveform.  Some  asymmetry  is  seen  in  the  wave  cycle  of 
the  peak-to-peak  ground  motion  for  the  events  on  November  23,  1984,  July  25,  1985,  and 
February  26,  1987,  where  the  first  half  of  the  cycle  is  of  slightly  lower  frequency  than  the 
second  half  of  the  cycle.  In  these  cases,  the  period  was  taken  as  twice  the  average  of  the 
cycle’s  half  periods.  The  explosion  on  April  3,  1987  was  clipped  at  NORESS,  and  the 
peak-to-peak  amplitude  was  estimated  by  graphically  extrapolating  the  waveforms  on 
each  side  of  the  clipped  cycle.  As  a  check,  an  mb  was  calculated  from  the  unclipped 
intermediate-period  channel  and  found  to  be  the  same  as  that  estimated  from  the  short- 
period  beam.  A  "P-factor"  of  3.33  (Veith  and  Claussen,  1972)  was  used  for  the  mb  calcu¬ 
lations.  Corrected  amplitudes,  periods,  calculated  magnitudes  and  residuals  are  summar¬ 
ized  in  Table  I. 

Figure  2  shows  a  plot  of  the  NORESS  mb  values  versus  the  NEIS  worldwide  average 
mb s,  as  well  as  the  residuals  (NORESS  mb  minus  NEIS  m4).  The  open  symbols  refer  to 
events  at  the  Degelen  test  site,  and  the  solid  symbols  refer  to  Shagan  tests.  Shagan  was 
the  site  of  the  seven  largest  tests  examined,  whereas  the  smaller  tests  were  conducted  at 
Degelen.  The  mean  residual  for  the  Degelen  tests  is  +0.21  and  the  mean  residual  for  the 
Shagan  tests  is  +1.00.  It  is  also  of  interest  to  notice  that  for  this  limited  dataset  the 
residuals  increase  with  increasing  magnitude  for  both. 


Table  I 


NEIS  Reported 

NORESS 

Date 

Lat. 

Long. 

7"  (sec) 

m6 

Residual 

Nov.  23,  1984 

49.897 

78.132 

El 

D 

8.2 

0.4 

4.64 

-0.06 

Dec.  16,  1984 

49.957 

78.862 

El 

S 

2806 

0.6 

7.00 

+0.90 

Apr.  25,  1985 

49.924 

78.969 

5.9 

s 

1582 

0.6 

6.75 

+0.85 

Jun.  15,  1985 

49.889 

78.881 

s 

3980 

0.6 

7.10 

+1.10 

Jun.  30,  1985 

49.861 

78.696 

6.0 

s 

2730 

7.07 

+1.07 

Jul.  20,  1985 

49.951 

78.829 

5.9 

s 

2492 

0.6 

6.91 

+1.01 

Jul.  25,  1985 

49.894 

78.150 

D 

47 

5.30 

+0.30 

Feb.  26,  1987 

49.787 

78.122 

5.2 

D 

89 

0.5 

5.58 

+0.38 

Mar.  12,  1987 

49.958 

78.834 

5.6 

s 

609 

0.6 

6.34 

+0.74 

Apr.  3,  1987 

50.59 

78.67 

6.1 

s 

1  1 

0.6 

7.3 

+1.3 

(1)  D  =  Degelen,  S  =  Shagan 

(2)  Peak-to-peak  amplitude  corrected  for  NORESS  instrument  response  at  period  T. 

(3)  Estimated  from  clipped  waveform.  See  text. 

Ringdal  (1986)  showed  that  there  is  a  bias  between  conventional  network  s  (as 
reported  by  NEIS)  and  maximum-likelihood  m^s,  especially  in  the  intermediate  magni¬ 
tude  range  (3.5  to  5.5).  Above  5.5,  the  two  magnitude  estimates  are  essentially  the  same. 
Thus,  the  NEIS  magnitudes  for  the  Shagan  explosions  studied  here  are  likely  unbiased, 
whereas  the  magnitudes  of  the  smaller  Degelen  shots  may  be  overestimated.  Using 
Ringdal’s  (1986)  results,  the  NEIS  magnitude  for  the  shot  on  Nov.  23,  1984,  should  be 
reduced  from  4.7  to  4.5,  giving  a  residual  for  NORESS  of  +0.14.  The  other  overestimated 
magnitude  would  be  that  of  the  July  25,  1985  shot,  which  would  be  reduced  from  5.0  to 
4.9,  yielding  a  NORESS  residual  of  +0.40.  The  mean  residual  for  the  Degelen  shots  thus 
becomes  +0.31  versus  the  previous  estimate  of  +0.21.  This  is,  however,  still  considerably 
lower  than  the  +1.00  mean  residual  for  the  Shagan  shots. 

5.2.3.  Conclusion 

Body  wave  magnitudes  were  calculated  for  ten  Eastern  Kazakh  explosions  using 
beams  of  NORESS  data.  The  NORESS  mb s  were  found  to  be  0.2  to  0.3  units  higher  than 
network  average  m^s  for  three  explosions  at  Degelen,  whereas  at  Shagan  the  bias  was 
+1.00.  It  remains  to  be  explained  why  the  difference  in  the  biases  should  be  so  large  for 
the  two  areas,  considering  that  they  are  only  50-70  km  apart.  It  should  also  be  recognized 
that  the  amplitude  versus  distance  curve  of  Veith  and  Claussen  (1972)  shows  a  large 
amount  of  scatter.  One  possibility  for  the  difference  in  residuals  may  be  the  focusing 
of  body  waves  by  Degelen  Mountain.  As  a  continuation  of  this  study,  more  Eastern 
Kazakh  events  will  be  studied  using  NORESS  data  and  the  results  combined  with  the 
magnitude  estimates  presented  here. 

Jay  J.  Pulli 
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Figure  1.  Beamed  P-waveforms  for  ten  Eastern  Kazakh  explosions  recorded  by  the  NORESS  array 
Body  wave  magnitudes  reported  by  NEIS.  The  waveform  for  event  1987:093  is  clipped. 
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Figure  2.  NORESS  m(s  versus  NEIS  m^s  as  well  as  residuals  for  the  ten  explosions 
studied. 
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5.3.  DETERMINATION  OF  AZIMUTHS  OF  SEISMIC  EVENTS  USING 
P  AND  RAYLEI^I  WAVE  POLARIZATION 


Well  determined  directions  of  arrival  of  P  or  other  seismic  phases  can  be  used  to 
help  obtain  an  initial  approximation  to  the  epicenter  of  an  event,  or  to  nelp  associate  an 
observed  phase  with  the  correct  epicenter  when  two  or  more  events  occur  at  about  the 
same  time.  The  objective  of  this  study  is  to  determine  if  recently  developed  analysis 
methods  can  improve  the  direction-finding  ability  of  three-component  seismic  stations, 
and  to  develop  automatic  methods  for  obtaining  this  information. 

A  polarization  filter  was  applied  to  a  data  set  extracted  from  the  GSETT  database 
consisting  of  waveforms  from  stations  RSNY,  RSON,  and  RSSD.  Three-component 
short-  and  long-period  data  were  obtained  from  67  and  56  events,  respectively,  with  a  typ¬ 
ical  magnitude  of  5.0.  Input  consisting  of  three-component  seismograms  (short-  and 
long-period)  are  first  band-pass  filtered  in  the  frequency  domain.  Next,  in  the  time- 
domain,  tapered  windows  are  applied  to  all  input  traces  where  the  time  resolution  is 
determined  by  the  window  length.  A  polarization  ellipsoid  is  computed  to  isolate  polar¬ 
ized  motion.  Parameters  of  the  ellipsoid  are  then  extracted  (Jurkevics,  1985). 

Table  I,  A  and  B  show  the  deviation  of  the  measured  azimuth  obtained  from  short- 
period  P- waves  and  long-period  Rayleigh  waves,  respectively,  from  the  calculated 
azimuth  to  the  epicenter.  As  noted  in  Table  I,  azimuths  measured  from  polarization 
analysis  of  both  P  and  Rayleigh  waves  are,  on  the  average,  good  estimates  of  the  direc¬ 
tion  to  the  epicenter.  On  the  other  hand,  individual  estimates  may  differ  significantly 
from  the  calculated  azimuth  of  the  epicenter.  For  P- waves,  RSSD  has  the  greatest 
dispersion  in  azimuth  of  arrival,  with  a  standard  deviation  of  79  *  .  RSON  has  a  much 
smaller  dispersion.  RSSD  shows  the  largest  deviation  for  Rayleigh  waves,  whereas  RSNY 
has  the  smallest.  Overall,  standard  deviations  for  all  three  stations  for  Rayleigh  waves 
are  quite  high.  Deviations  may  be  attributed  to  misorientation  of  the  instruments;  failure 
of  the  theoretical  wave  model  of  the  study,  or  may  be  due  to  inhomogenities  of  the  earth. 

Additional  statistical  analysis  is  needed  to  be  more  certain  of  these  preliminary 
results.  This  may  include  dependence  of  accuracy  on  the  signal-to-noise  ratio,  and  the  rec- 
tilinearity  or  planarity  of  the  waveform. 


Antoinette  Campanella 
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Table  I.  Azimuthal  Anomalies  of  Selected  Events  for  RSNY,  RSON  and  RSSD 


(A  -  P  waves) 


Short-Period 

Arithmetic 
Mean  (deg.) 

Standard 
Deviation  (deg.) 

Median  (deg.) 

All  Stations 
(absolute  values) 

10.31 

52.56 

11.43 

RSNY 

-16.81 

30.61 

16.51 

RSON 

2.03 

21.78 

0.10 

RSSD 

12.10 

79.29 

17.67 

(B  -  Rayleigh  waves) 


Long-Period 

Arithmetic 
Mean  (deg.) 

Standard 
Deviation  (deg.) 

Median  (deg.) 

All  Stations 
(absolute  values) 

7.19 

62.25 

7.80 

RSNY 

-7.76 

57.66 

10.19 

RSON 

6.63 

57.56 

2.08 

RSSD 

7.19 

69.83 

11.14 
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5.4.  INVENTORY  OF  RUSSIAN  FILM  DATA  AT  THE  CENTER 


In  a  preliminary  study  of  seismic  signals  recorded  in  the  U.S.S.R.,  Ryall  (1986)  digi¬ 
tized  tracings  of  P- waves  from  copies  of  Russian  short-period  recordings  of  several  earth¬ 
quakes  in  southeastern  California.  Analysis  of  the  Soviet  data,  and  comparison  of  magni¬ 
tude  residuals  with  those  obtained  in  previous  studies,  led  to  an  estimate  of  magnitude 
bias  between  the  Nevada  Test  Site  and  a  number  of  station  sites  in  the  U.S.S.R.,  includ¬ 
ing  a  station  (SEM)  at  Semipalatinsk,  about  100  km  from  the  Eastern  Kazakh  test  site. 

In  addition  to  the  original  data  sample  we  now  have  copies  on  hand  of  long-  and 
short-period,  three-component  recordings  from  about  40  Soviet  stations,  either  requested 
by  other  researchers  through  World  Data  Center  A  or  supplied  as  part  of  the  Interna¬ 
tional  Data  Exchange.  This  data  set  includes  film  copies  of  recordings  of  earthquakes  in 
China,  Japan,  the  Sea  of  Okhotsk,  Alaska,  Austria,  Iran,  Kamchatka,  the  Caucasus, 
Romania,  and  regions  farther  from  the  U.S.S.R.  ( Figure  1).  Typically,  data  are  available 
for  a  dozen  or  more  stations  for  a  given  earthquake,  and  the  quality  of  the  recordings 
varies  from  very  poor  to  very  good.  Station  calibration  constants  are  either  marked  on 
the  records  or  available  from  the  Soviet  bulletins  or  other  published  reports. 

These  records  contain  direct  information  that  can  be  used  to  estimate  characteristics 
of  body  and  surface  wave  propagation  in  the  U.S.S.R.  Studies  which  may  be  carried  out 
with  these  data  include  the  determination  of  receiver  structure  under  each  station  from 
three-component  data,  magnitude  bias  versus  azimuth,  P-  and  5-wave  residuals  versus 
azimuth,  and  a  comparison  of  these  characteristics  with  the  tectonics  of  the  U.S.S.R. 
Crust  and  upper-mantle  structure  as  well  as  attenuation  can  be  estimated  using  the  two- 
station  transfer  function  method  applied  to  surface  wave  paths  crossing  the  area.  If  local 
and  regional  events  happen  to  be  present  on  the  records,  the  data  can  be  used  to  estimate 
local  attenuation  and  site  effects  at  the  station.  The  results  of  such  analyses  can  be  used 
to  increase  our  confidence  in  network  evaluations  by  programs  such  as  SNAP-D,  and 
should  make  a  significant  contribution  to  resolving  controversies  stemming  from  differ¬ 
ences  in  currently  assumed  propagation  characteristics  within  the  U.S.S.R. 

As  a  fi-st  step  toward  creation  of  a  database  containing  a  representative  sample  of 
signals  digitized  from  this  record  collection,  we  have  compiled  an  inventory  of  the  data  on 
hand.  Table  I  lists  50  events  for  which  an  adequate  amount  of  good-quality  data  is  avail¬ 
able  for  analysis,  and  Table  II  lists  28  seismic  stations  which  supplied  most  of  the  data. 
Figure  2  shows  the  locations  of  these  stations  on  a  polar  projection  of  the  U.S.S.R.  cen¬ 
tered  on  latitude  66  *  N,  longitude  90  *  E.  Table  III  indicates  the  quality  of  short-period 
recordings  available  for  the  50  events,  with  solid  circles  indicating  "good"  to  "very  good" 
quality,  open  circles  "fair"  and  x’s  "poor"  quality.  Seismograms  classified  as  "very  good" 
had  traces  on  which  P  and  other  phases  were  clearly  recorded.  Those  rated  as  "good" 
had  P~ wave  onsets  that  were  clear  for  at  least  several  seconds,  but  other  parts  of  the 
record  may  have  been  unusable.  "Fair"  recordings  typically  had  P  onsets  that  were 
fogged  or  dim  but  probably  traceable,  but  other  parts  of  the  signal  may  be  poor.  "Poor" 
signals  were  generally  unusable  as  the  result  of  under-  or  overexposure,  intermixed  traces, 
smudging,  fogging,  etc.  Of  a  total  of  457  sets  of  records,  41%  were  rated  as  "poor,"  23% 
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III.  Circle  has  radius  of  90  * ,  is  centered  on 


Table  I.  Events  for  which  film  data  are  available. 


Date 

Time 

Lat. 

Long. 

650925 

1547:56 

41.52N 

75.04E 

670211 

0927:34 

52.16N 

106.48E 

691124 

1723:19 

37.16N 

71.64E 

700828 

0745:59 

6.26S 

154.62E 

700828 

0944:57 

52.27N 

105.90E 

700830 

1746:08 

52.36N 

151.64E 

700905 

0752:27 

52.28N 

151.49E 

710129 

2158:03 

51.69N 

150.97E 

710328 

2047:16 

41.42N 

79.20E 

710824 

0226:14 

05.28S 

151.54E 

710403 

0449:03 

32.16N 

94.99E 

710510 

1451:45 

42.85N 

71.29E 

710615 

2204:13 

41.39N 

79.18E 

740510 

1925:17 

28.20N 

104.00E 

740511 

0056:56 

1.89N 

126.47E 

741003 

1111:56 

36.50N 

70.83E 

741003 

1421:29 

12.24S 

77.58W 

741008 

0950:58 

17.37N 

61.99W 

750526 

0911:51 

35.98N 

17.56W 

751011 

1435:15 

24.91S 

175.18W 

751226 

1556:39 

16.22S 

172.47W 

760114 

1556:33 

29.48S 

177.57W 

760204 

0930:35 

14.90N 

90.53W 

760506 

2000:12 

46.35N 

1S.26E 

760727 

1942:54 

39.56N 

117.87E 

760728 

0058:44 

39.40N 

117.76E 

760728 

0948:23 

39.44N 

118.03E 

760728 

1045:37 

39.71N 

118.37E 

760816 

1611:05 

6.22N 

124. 10E 

761124 

1222:16 

39.05N 

44.04E 

770304 

1921:54 

45.83N 

26.72E 

770402 

0715:22 

16.79S 

172.02W 

770404 

1906:35 

02.78S 

102. 26E 

770420 

2342:53 

09.94S 

160.48E 

770421 

0424:09 

10.00S 

160.77E 

770822 

1208:33 

22.91S 

175.74W 

770819 

0608:54 

11.16S 

118. 41E 

790228 

2127:06 

60.74N 

141.55W 

791127 

1710:33 

34.08N 

59.79E 

791212 

0759:04 

1.62N 

79.34W 

800525 

1633.44 

37.60N 

118.80W 

800527 

1451:00 

37.49N 

118.78W 

800708 

2319:19 

12.49S 

166. 37E 

800717 

1942:23 

12.48S 

166.06E 

801010 

1225:22 

36.16N 

1.40E 

801108 

1027:32 

41.15N 

124. SOW 

810901 

0929:31 

15.08S 

173.12W 

830318 

0905:50 

4.89S 

153.57E 

831130 

1746:00 

6.85S 

72.04E 

840207 

2133:20 

9.96S 

160. 49E 

Depth  m4  Remarks 


12 

5.2 

Hindu  Kush  &  Pamir 

26 

5.3 

Alma  Ata  to  Lake  Baikal 

113 

5.8 

Pamir-Hindu  Kush 

5.8 

Solomon  Islands 

30 

4.9 

Lake  Baikal  Region 

643 

6.5 

Sea  of  Okhotsk 

5.7 

Sea  of  Okhotsk 

515 

6.0 

Sea  of  Okhotsk 

14 

(5.8) 

Kirgiz-Sinkiang  Border 

69 

(5.7) 

New  Britain  Region 

27 

(5.6) 

Tibet 

14 

5.6 

Kirgiz 

96 

(5.3) 

Kirgiz-Sinkiang  Border 

17 

5.8 

Szechwan  Province 

50 

5.4 

Molucca  Passage 

186 

4.9 

Hindu  Kush  Region 

9 

6.2 

Near  Coast  of  Peru 

41 

8.4 

Leeward  Islands 

34 

6.5 

N.  Madeira 

11 

6.4 

S.  Tonga 

33 

6.3 

Samoa 

29 

6.4 

E.  Kermadec  Islands 

44 

5.4 

Guatemala 

12 

5.9 

Austria  -  N.E.  Italy 

10 

6.1 

N.E.  China 

12 

4.9 

N.E.  China 

35 

4.9 

N.E.  China 

37 

6.1 

N.E.  China 

8 

6.4 

Phillipines- Mindanao 

10 

6.1 

N.W.  Iran-USSR  Border 

86 

6.1 

Romania 

33 

6.4 

N.E.  Tonga 

126 

5.3 

S.  Sumatra 

33 

6.1 

Solomon  Islands 

31 

6.4 

Solomon  Islands 

69 

6.3 

Tonga 

78 

6.8 

S.  of  Sumbawa 

22 

6.2 

S.E.  Alaska 

9 

6.2 

Iran 

28 

6.2 

Near  Coast  of  Ecuador 

? 

6.1 

Mammoth  Lakes,  Calif. 

? 

5.7 

Mammoth  Lakes,  Calif. 

29 

5.9 

Santa  Cruz  Islands 

29 

5.6 

Santa  Cruz  Islands 

0 

6.3 

Algeria 

5 

5.9 

Near  Coast  of  N.  Calif. 

20 

6.5 

Tonga 

94 

6.5 

New  Ireland  region 

10 

6.5 

Chagos  Archipelago 

4 

6.5 

Solomon  Islands 

Table  II.  Station  Information. 


Name 

Code 

Lat.  N. 

Long.  E. 

Andizhan 

AND 

40.75 

72.37 

Apatity 

APA 

67.55 

33.33 

Bakuriani 

BKR 

41.73 

43.52 

Boday bo 

BDB 

57.85 

114.18 

Dushanbe 

DSH 

38.57 

68.77 

El’tsovka 

ELT 

53.25 

86.27 

Frunze 

FRU 

42.83 

74.62 

Goris 

GRS 

39.50 

46.33 

Irkutsk 

IRK 

52.27 

104.32 

Iul’tin 

ILT 

67.83 

178.80 

Kheys 

KHE 

80.62 

58.05 

Kirovabad 

KRV 

40.65 

46.33 

Magadan 

MAG 

59.55 

150.80 

Makhachkala 

MAK 

43.02 

47.43 

Mondy 

MON 

51.68 

100.98 

Murgab 

MUR 

38.37 

73.93 

Novosibirsk 

NVS 

55.03 

82.92 

Obninsk 

OBN 

55.17 

36.45 

Petropavlovsk 

PET 

53.02 

158.65 

Przheval’sk 

PRZ 

42.48 

78.40 

Semipalatinsk 

SEM 

50.40 

80.25 

Simferopol’ 

SIM 

44.95 

34.12 

Sochi 

soc 

43.58 

39.72 

Talgar 

TLG 

43.23 

77.23 

Tiksi 

TIK 

71.63 

128.87 

Uzhgorod 

UZH 

48.63 

22.30 

Y  uzhno-Sakhalinsk 

YSS 

47.02 

142.72 

Zakamensk 

ZAK 

50.38 

103.28 
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"fair"  and  36%  "good"  to  "very  good."  Long-period  seismograms  were  generally  clearer 
than  the  short-period  records  characterized  by  Table  III,  but  an  inventory  of  long-period 
data  has  not  been  made. 

Station  SEM,  because  of  its  proximity  to  the  Eastern  Kazakh  test  range,  is  of  partic¬ 
ular  interest  to  questions  of  propagation  bias.  That  station  had  good  quality  records  for 
seven  of  the  50  events,  fair  for  another  three,  and  poor  for  nine.  In  addition,  SEM  had 
good  recordings  of  two  events  in  1965  for  which  the  ISC  Bulletin  did  not  list  locations. 

The  Center  plans  to  begin  hand-digitizing  the  best  of  the  Soviet  short-period  records 
over  the  next  few  months,  and  installing  the  digitized  traces  in  a  special  database.  Pro¬ 
gress  on  this  effort  will  be  reported  in  a  Center  Newsletter,  which  will  be  issued  quarterly 
starting  in  June. 


Alan  S.  Ryall,  Jr. 
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5.5.  STATUS  OF  NORESS  AND  OTHER  WAVEFORM  DATABASES 


5.5.1.  Introduction 

A  major  goal  of  the  Center  program  has  been  to  facilitate  the  work  of  other  contrac¬ 
tors  in  key  areas  of  the  nuclear  monitoring  research  program,  with  the  Center  supplying 
data  and  computer  services  to  external  users.  The  Center  currently  maintains  or  is 
preparing  a  sizeable  number  of  databases.  Some  of  these  represent  a  complete  set  of  infor¬ 
mation  of  some  type  and  the  contents  are  fixed,  while  others  are  undergoing  development 
on  a  continual  basis  as  new  data  of  the  same  type  arrive. 

5.5.2.  Status  of  Prrnary  Databases 

All  Center  databases  are  on  HUGO,  and  are  accessed  through  the  INGRES  rela¬ 
tional  database  management  system.  HUGO  currently  allocates  1,084  Mbytes  of  disk 
space  for  databases,  of  which  945  Mbytes  are  filled.  Status  of  the  primary  databases  are 
given  in  Table  I.  Note  in  relation  to  the  waveform  databases  that  the  term  "complete” 
means  that  the  Center  has  all  data  available  for  a  particular  source;  "incomplete"  means 
that  additional  data  may  be  available  from  another  agency.  Some  of  the  waveform  data 
at  the  Center  is  in  the  form  of  waveform  segments  for  particular  events  (GDSN, 
WAKE,  SAFRICA,  HDWWSSN,  NORSAR),  while  the  remaining  data  sets  con¬ 
tain  both  continuous  and  segmented  data. 


Table  I.  Status  of  Primary  Center  Databases 


Database 

Coverage 

No.  Tapes 

Mbytes  (disk) 

Parameter  Data 

EVENTS 

1177  BC  to  current 

0 

52.544 

ISC 

Nov.  1,  1978  to  Dec.  SI,  1978 

0 

19.760 

ARRIVALS 

May  16,  1982  to  current 

0 

175.196 

GSETT 

Oct.  15,  1984  to  Dec.  14,  1984 

0 

21.524 

Waveforms-complete 

GDSN 

Jan.  1,  1976  to  Dec.  SI,  1986 

104  S 

256.166 

RSTN 

Jan.  1,  1982  to  current 

2975 

153.024 

(dp/pdp  Jan.  15,  1984  to  current) 

0 

59.640 

WAKE 

Sep.  8,  1982  to  Dec.  10,  1986 

138 

3.916 

SAFRICA 

Mar.  15,  1986  to  Oct.  11,  1986 

1 

0.448 

NORESS  (on-line) 

Oct.  29,  1986  to  current 

384 

20.800 

NORESSHF 

Dec.  1,  1985  to  Dec.  31,  1986 

183 

1.360 

Waveforms- incomplete 

NORESS 

Oct.  25,  1984  to  Oct.  28,  1986 

172 

19.116 

HDWWSSN 

Sep.  IS,  1963  to  Dec.  27,  1981 

4 

1.752 

NORSAR 

Nov.  6,  1976  to  Dec.  IS,  1986 

100 

18.920 

Additional  sets  of  waveforms  that  are  being  processed  for  inclusion  as  official  Center  data¬ 
bases  include  Phase  I  and  Phase  II  data  from  stations  being  installed  and  operated  in 
Kazakhstan,  U.S.S.R.,  in  a  cooperative  effort  by  the  Natural  Resources  Defense  Council 
and  the  Soviet  Academy  of  Sciences,  and  a  large  number  of  recordings  of  NTS  explosions 
made  under  the  Long  Range  Seismic  Measurements  program  and  the  Special  Data  Collec¬ 
tion  Systems  program.  Work  is  also  progressing  on  compilation  of  classified  and  unclassi¬ 
fied  versions  of  an  EXPLOSION  database,  and  on  the  archiving  of  instrument  responses 
for  all  waveform  data  available  at  the  Center. 

5.5.3.  NORESS  Research  Data  Sets 

With  the  DARPA/GSD  program  now  emphasizing  advanced  research  utilizing  data 
from  small  seismic  arrays,  the  Center  has  devoted  considerable  effort  to  compilation  of  a 
representative  set  of  NORESS  waveform  segments  for  local  and  regional  events  over  a 
range  of  azimuths  from  the  NORESS  array,  plus  waveforms  for  all  underground  nuclear 
explosions  recorded  by  the  array.  These  data  are  archived  in  three  separate  data  sets, 
described  below  and  summarized  in  Tables  II-IV.  Center  users  may  request  copies  of 
these  data  sets  on  tar  tapes,  by  contacting  Richard  Baumstark. 

The  Center  research  staff  has  selected  recordings  of  88  small  earthquakes  and  mine 
blasts  in  Scandinavia  and  the  western  U.S.S.R.  as  representing  the  types  of  events 
recorded  over  a  range  of  distance  out  to  15  *  and  a  range  of  azimuth  around  the  array. 
Epicentral  information  for  this  set  of  waveforms,  which  are  currently  retained  on  disk  on 
one  of  the  Sun  computers  (maut),  is  listed  in  Table  II  together  with  event  type  taken  from 
local  station  bulletins.  As  the  table  indicates,  the  set  includes  34  earthquakes,  38  mine 
blasts,  and  16  events  for  which  event  type  was  not  specified  in  the  local  bulletins.  Ten  of 
the  latter  occurred  off  the  southwest  coast  of  Norway  on  November  20-21,  1985,  and  were 
probably  earthquakes.  Locations  of  the  88  events  are  shown  at  the  top  of  Figure  1 ,  with 
the  circle  drawn  at  a  distance  of  15  *  from  NORESS. 


Table  II.  Events  Selected  by  Center  Research  Stall 


Name 

Mo-Dy 

HrMn:Sec 

Lat 

Lon 

Dist 

Depth 

mb 

Area 

Type 

85213 

08-01 

11:17:35.8 

45.82 

26.65 

17.39 

122 

4.7 

Romania 

E 

85298 

10-25 

12:03:47. 

59.30 

28.10 

8.42 

0 

2.3 

SW  of  Leningrad 

X 

85300 

10-27 

04:38:43 

61.12 

4.92 

3.26 

0 

2.8 

W  Coast  Norway 

E 

85312 

11-08 

14:18:54.6 

58.34 

6.43 

3.54 

0 

2.4 

SW  Norway 

X 

85313 

11-09 

14:42:46. 

57.80 

7.20 

3.69 

0 

2.1 

SW  Norway 

X 

85313b 

11-09 

18:20:48. 

62.00 

7.70 

2.24 

0 

2.0 

Norway  due  NT 

T 

85317 

11-13 

16:32:10. 

58.30 

6.40 

3.54 

0 

1.8 

SW  Norway 

X 

85317b 

11-13 

12:07:48. 

59.30 

28.10 

8.42 

0 

2.3 

SW  of  Lenin 

X 

85324 

11-20 

22:10:44.2 

57.61 

5.67 

4.35 

0 

2.3 

SW  Norway 

T 

85324b 

11-20 

22:24:38.1 

57.66 

5.72 

4.30 

0 

2.2 

SW  Norway 

? 

85324c 

11-20 

22:57:10.0 

57.64 

5.62 

4.35 

0 

2.3 

SW  Norway 

T 

Table  IL  Event*  Selected  by  Center  Bee* arch  Staff  —  (Continued) 


Nam* 

Mo-Djr 

Hr.Mn:S«c 

Let 

Loa 

Diet 

Depth 

mb 

Area 

Type 

85324d 

11-20 

23:10:47.5 

57.88 

5.35 

4.43 

0 

2.3 

SW  Norway 

? 

85324* 

11-20 

23:17:28.9 

57.89 

5.49 

4.38 

0 

2.3 

SW  Norway 

? 

8S324f 

11-20 

23:23:10.0 

57.50 

5.62 

4.45 

0 

2.2 

SW  Norway 

? 

85324g 

11-20 

23:28:23.1 

57.58 

5.49 

4.44 

0 

2.2 

SW  Norway 

? 

85325 

11-21 

14:18:13. 

59.80 

8.20 

1.91 

0 

1.4 

SW  Norway 

? 

85325b 

11-21 

14:48:07. 

54.80 

8.50 

8.53 

0 

23 

SW  Norway 

? 

85325c 

11-21 

09:18:30. 

58.37 

12.38 

2.41 

-.0 

SE  Norway? 

T 

85327 

11-23 

13:08:18. 

59.50 

25.00 

8.84 

0 

2.1 

NW  Estonia 

X 

85331 

11-27 

04:53:32.8 

59.73 

05.71 

3.08 

15 

2.8 

W  Coaet  Norway 

E 

85344 

12-10 

12:05:39. 

59.40 

28.50 

8.58 

0 

2.2 

SW  of  1»— in  grad 

X 

85357 

12-23 

02:35:08.3 

80.38 

01.90 

4.77 

15 

2.3 

Shetland  b 

E 

85358 

12-24 

12:37:57. 

59.80 

22.50 

5.53 

1.9 

S  tip  of  Finland 

? 

SS3S9 

12-25 

13:19:01. 

58.70 

18.00 

7.59 

0 

2.8 

Estonia 

? 

85381 

12-27 

12:18:08 

59.40 

28.50 

8.58 

0 

2.4 

SW  of  Leningrad 

X 

85385 

12-31 

08:57:17. 

73.31 

6.82 

12.77 

15 

4.8 

Cropland  S*a 

E 

66003 

01-03 

14:58:41. 

81.90 

30.80 

9.23 

0 

2.5 

NW  of  Leningrad 

X 

88007 

01-07 

14:14:28. 

58.34 

6.43 

3.54 

0 

2.2 

SW  Norway 

X 

88009 

01-09 

08:18:43. 

54.70 

19.50 

7.40 

2.7 

Poland 

? 

88017 

01-17 

14:11:01. 

84.34 

6.43 

3.54 

0 

2.3 

SW  Norway 

X 

88019 

01-19 

04:59:21. 

85.00 

12.13 

4.29 

3.0 

N  Sweden 

E 

88020 

01-20 

23:38:28. 

50.19 

12.37 

10.58 

4.9 

Germany 

E 

88021 

01-21 

08:55:40. 

55.30 

13.80 

5.58 

2.5 

S  tip  of  Sweden? 

? 

88031 

01-31 

12:10:15. 

59.30 

28.10 

8.43 

0 

3.2 

SW  Leningrad 

X 

88035 

02-04 

12:14:59. 

59.50 

28.50 

7.58 

0 

2.8 

NX  Female 

X 

88035b 

02-04 

12:58:59. 

59.4 

24.80 

9.97 

0 

2.5 

NW  Batenin 

X 

1603  fie 

02-04 

14:22:57. 

•9.30 

24.40 

9.91 

0 

2.8 

NW  Eetetua 

X 

66036 

02-08 

17:53:19. 

92.81 

4.88 

3.80 

4.7 

Norwegian  3c* 

E 

88037 

02-08 

18:39:55. 

87.10 

10.80 

7.53 

0 

2.7 

N  Swtdm 

X 

88037b 

02-08 

12:22:04. 

59.30 

28.10 

8.42 

0 

2.7 

SW  Leningrad 

X 

88038b 

02-07 

11:00:01. 

84.70 

30.70 

9.84 

0 

3.1 

EareSaa  MB 

X 

88041 

02-10 

12:41:48. 

59.40 

28.50 

8.58 

0 

3.5 

SW  of  Lenta 

X 

66Q  H 

02-14 

14:13:19. 

58.34 

9.43 

3.54 

0 

2.4 

SW  Norway 

X 

88045b 

02-14 

17:54:04. 

58.34 

9.43 

3.54 

0 

2.3 

SW  Norway 

X 

88045c 

02-14 

12:10:21. 

59.40 

28.50 

8.58 

0 

1.7 

SW  of  Leningrad 

X 

880458 

02-14 

18:44:08. 

97.10 

20.80 

7.53 

0 

2.9 

N  Sweden 

X 

88049 

02-18 

10:48:19. 

59.30 

27  20 

8.00 

0 

2.8 

NX  Estonia 

X 

88049b 

02-18 

12:45:50. 

94.70 

30.70 

•  84 

0 

2.8 

Karelian  SSR 

X 

88087 

02-28 

02:11:44. 

92.79 

05  29 

3.80 

15 

2.5 

W  Coaet  Norway 

X 

88082 

03-03 

07:28:08 

43  70 

31  40 

20  81 

45 

4.4 

Stack  Sea 

E 

66064 

03-08 

14:18:31. 

99.3 

31  7 

7.19 

1 

-  0 

N  Sweden* 

? 

660646 

03-08 

12:13:19. 

59  50 

28  50 

758 

0 

2.9 

NX  Estonia 

X 

88084c 

03-08 

13.-02:05. 

80  83 

02.59 

4  41 

15 

2  1 

W  Coaet  Norway 

X 

1-14 


Table  II.  Event*  Selected  by  Center  Research  Staff  —  (Continued) 


Name 

Mo-Dy 

Hr.Mn:Sec 

Lat 

Lon 

Dist 

Depth 

mb 

Area 

Type 

86067 

03-08 

16:21:17. 

61.67 

02.58 

4.44 

IS 

2.4 

W  Coast  Norway 

E 

86069 

03-10 

12:02:09. 

59.30 

28.10 

8.42 

0 

2.6 

SW  of  Leningrad 

X 

86069b 

03-10 

04:20:04. 

62.81 

04.91 

3.78 

IS 

2.5 

W  Coast  Norway 

E 

86070 

03-11 

12:02:28. 

59.30 

28.10 

8.42 

0 

2.6 

SW  of  Leningrad 

X 

86071 

03-12 

11:07:21. 

59.50 

26.50 

7.57 

0 

2.5 

NE  Estonia 

X 

86071b 

03-12 

12:01:38. 

59.40 

28.50 

8.58 

0 

2.5 

SW  of  Leningrad 

X 

86078 

03-19 

12:06:40. 

59.40 

28.50 

8.58 

0 

2.6 

SW  of  Leningrad 

X 

86089 

03-30 

03:22:37. 

61.66 

04.53 

3.52 

15 

2.2 

W  Coast  Norway 

E 

86091 

04-01 

09:56:53. 

56.42 

12.10 

4.35 

IS 

3.6 

Sweden 

E 

86094 

04-04 

22:42:30. 

71.08 

08.35 

10.47 

IS 

4.6 

Norwegian  Sea 

E 

86097 

04-07 

00:34:37. 

61.84 

04.88 

3.40 

15 

2.3 

W  Coast  Norway 

E 

86108 

04-18 

00:44:13. 

59.22 

01.42 

5.30 

15 

2.4 

Shetland  Is 

E 

861S4 

06-03 

14:30:04. 

61.46 

04.08 

3.69 

IS 

2.8 

W  Coast  Norway 

E 

86158 

06-04 

09:06:31 

61.50 

30.40 

9.16 

0 

3.3 

N  of  Leningrad 

X 

86163 

06-12 

09:30:55. 

61.50 

30.40 

9.16 

0 

3.1 

N  of  Leningrad 

X 

86166 

06-15 

15:01:07. 

61.67 

03.85 

3.84 

IS 

3.0 

W  Coast  Norway 

E 

86168 

06-17 

12:12:07. 

59.4 

28.5 

8.58 

0 

2.6 

SW  of  Leningrad 

X 

86169 

06-18 

11:05:08. 

59.4 

28.5 

8.58 

0 

2.5 

SW  of  Leningrad 

X 

86170 

06-19 

03:55:08. 

59.31 

06.54 

2.89 

0 

2.4 

W  Coast  Norway 

E 

86171 

06-20 

22:07:53. 

61.47 

03.92 

3.77 

15 

2.0 

W  Coast  Norway 

E 

86177 

06-26 

04:06:21. 

61.88 

05.10 

3.31 

15 

2.4 

W  Coast  Norway 

E 

86178 

06-27 

03:49:46. 

59.28 

06.76 

2.81 

15 

2.5 

W  Coast  Norway 

E 

86185 

07-04 

11:13:27. 

59.3 

28.1 

8.42 

0 

2.6 

SW  of  Leningrad 

X 

86195 

07-14 

13:50:32. 

58.35 

13.82 

2.66 

10 

4.0 

S  Sweden 

E 

86195b 

07-14 

14:30:27. 

61.10 

29.90 

8.95 

0 

2.9 

N  Finland 

X 

86195c 

07-14 

15:02:19. 

69.30 

34.40 

12.83 

0 

2.9 

Kola  Peninsula 

X 

86223 

08-10 

08:01:04. 

59.99 

05.34 

3.17 

IS 

1.7 

W  Coast  Norway 

E 

86228 

08-16 

04:24:36. 

62.82 

04.98 

3.75 

IS 

2.5 

W  Coast  Norway 

E 

86244 

09-01 

22:11:26. 

60.82 

02.93 

4.22 

15 

3.5 

W  Coast  Norway 

E 

86273 

09-30 

20*2:47. 

80.79 

04.23 

3.59 

15 

2.4 

W  Coast  Norway 

E 

86283 

10-10 

19:56:31. 

61.97 

02.33 

4.60 

15 

2.3 

W  Coast  Norway 

E 

86299 

10-26 

11:45:06. 

61.46 

03.29 

4.07 

15 

2.6 

W  Coast  Norway 

E 

86299b 

10-28 

11:57:03. 

61.72 

03.27 

4.12 

IS 

2.6 

W  Coast  Norway 

E 

88302 

10-29 

21:05:01. 

60.81 

03.04 

4.17 

IS 

2.4 

W  Coast  Norway 

E 

86327 

11-23 

03:30:32. 

73.74 

09.08 

13.10 

10 

4.7 

Green  land  Sea 

E 

86346 

12-12 

16:33:30. 

72.96 

04.80 

12.55 

10 

4.7 

Norwegian  Sea 

E 

(•II 


Table  111.  Events  Selected  by  S  AIC,  La  Jolla  Group 


Lat 

Lon 

Date 

Start  Time 

Depth 

mb 

59.30 

6.95 

Nov.  06,  1985 

14:50:50.2 

10 

2.3 

58.00 

6.60 

Nov.  13,  1985 

14:11:05.8 

0 

1.9 

59.70 

7.70 

Nov.  21,  1985 

13:16:47.2 

0 

1.7 

58.80 

4.80 

Nov.  21,  1985 

14:06:44.4 

0 

2.0 

56.70 

7.30 

Nov.  21,  1985 

15:04:34.7 

0 

2.3 

58.30 

4.80 

Nov.  21,  1985 

15:47:54.1 

0 

2.1 

61.77 

4.49 

Nov.  30,  1985 

19:05:09.4 

10 

3.0 

60.60 

10.10 

Dec.  04,  1985 

14:00:01.2 

0 

1.2 

60.20 

5.34 

Dec.  07,  1985 

14:15:32.3 

10 

0.0 

58.80 

5.70 

Dec.  07,  1985 

14:39:09.0 

0 

2.0 

58.30 

6.10 

Dec.  31,  1985 

13:36:49.6 

0 

2.1 

58.20 

6.90 

Jan.  31, 1986 

14:17:35.7 

0 

1.9 

58.00 

9.60 

Feb.  03,  1986 

15:52:59.5 

0 

2.8 

62.40 

5.80 

Feb.  05,  1986 

23:35:53.9 

0 

2.1 

62.30 

5.70 

Feb.  06,  1986 

06:20:05.4 

0 

1.8 

62.10 

6.50 

Mar.  25,  1986 

09:05:33.9 

0 

1.8 

58.70 

5.80 

Apr.  04,  1986 

13:12:43.9 

0 

1.9 

60.30 

4.90 

Apr.  08,  1986 

11:27:36.1 

0 

1.7 

58.40 

10.55 

Apr.  09,  1986 

08:16:24.9 

10 

2.5 

58.40 

10.60 

Apr.  09,  1986 

08:27:01.7 

0 

2.4 

60.30 

5.39 

Apr.  16,  1986 

11:51:11.5 

10 

1.8 

57.90 

7.40 

Apr.  16,  1986 

13:14:36.8 

0 

2.1 

60.80 

4.20 

Apr.  28,  1986 

15:52:49.9 

0 

2.5 

58.30 

4.80 

Apr.  30,  1986 

06:04:47.7 

0 

1.9 

56.20 

1.30 

Apr.  30,  1986 

06:21:05.7 

0 

2.7 

59.31 

6.88 

Apr.  30,  1986 

10:18:47.3 

10 

2.3 

60.36 

5.08 

May  02,  1986 

10:44:57.0 

10 

1.8 

59.70 

10.90 

May  05,  1986 

15:38:10.8 

0 

0.7 

63.30 

6.40 

May  17,  1986 

16:00:40.5 

0 

2.4 

59.70 

7.20 

May  27,  1986 

18:36:14.0 

0 

2.1 

59.30 

7.30 

May  28,  1986 

17:51:56.8 

0 

2.4 

59.60 

4.90 

Jun. 03, 1986 

11:03:23.9 

0 

2.1 

58.50 

6.50 

Jun. 06, 1986 

13:14:27.9 

0 

1.7 

58.90 

5.90 

Jun. 23, 1986 

13:12:54.0 

0 

1.8 

59.30 

7.30 

Jul.  10,  1986 

20:10:28.2 

0 

2.3 

62.50 

7.60 

Jul.  12,  1986 

13:37:46.6 

0 

1.9 

57.70 

14.30 

Jul.  14,  1986 

14:45:18.0 

0 

3.5 

66.80 

6.80 

Jul.  15,  1986 

18:45:32.9 

0 

3.5 

59.23 

7.23 

Jul.  16,  1986 

17:49:24.8 

10 

2.6 

5-20 


Table  HI 

.  Events  Selected  by  SAIC, 

La  Jolla  Group 

—  (Continued) 

Lat 

Lon 

Date 

Start  Time 

Depth 

mb 

59.30 

7.40 

Jul.  23,  1986 

20:47:09.8 

0 

2.2 

59.30 

7.30 

Jul.  29,  1986 

13:13:40.5 

0 

2.3 

59.50 

7.00 

Jul.  30,  1986 

17:59:39.0 

0 

2.4 

60.00 

10.60 

Jul.  30,  1986 

22:30:13.5 

0 

2.1 

58.30 

6.70 

Aug.  14,  1986 

13:14:39.1 

0 

1.9 

59.30 

7.40 

Aug.  14,  1986 

14:39:57.3 

0 

2.4 

59.50 

7.00 

Sep.  02, 1986 

12:53:50.8 

0 

2.1 

59.20 

7.40 

Sep.  09, 1986 

17:55:57.5 

0 

2.4 

59.20 

16.00 

Sep.  20,  1986 

22:14:59.3 

0 

3.5 

58.50 

6.40 

Oct.  01,  1986 

14:15:09.6 

0 

1.9 

57.80 

7.20 

Oct.  09,  1986 

14:13:51.5 

0 

2.0 

61.30 

6.50 

Oct.  12,  1986 

14:53:15.0 

0 

1.5 

59.60 

10.80 

Oct.  30,  1986 

09:21:15.5 

0 

1.2 

62.60 

6.30 

Nov.  01,  1986 

14:55:03.8 

0 

2.4 

58.90 

12.70 

Nov.  02,  1986 

07:48:02.5 

0 

3.5 

57.70 

8.70 

Nov.  13,  1986 

08:00:29.7 

0 

1.8 

I 


6*37 


Table  IV.  Nuclear  Explosions  Recorded  at  NORESS 


I 
1 

I 

LAT  LON  DATE  TIME  MB  GRN 


1  73.37 

54.96 

Oct.  25,  1984 

06:29:57.7 

5.9 

648 

49.95 

78.83 

Oct.  27,  1984 

01:50:10.6 

6.2 

329 

■  46.86 

48.10 

Oct.  27,  1984 

05:59:57.1 

5.0 

357 

0  46.84 

48.08 

Oct.  27,  1984 

06:04:56.7 

5.0 

357 

37.00 

-116.02 

Nov.  10,  1984 

16:40:00.0 

4.5 

41 

n  49.90 

78.13 

Nov.  23,  1984 

03:55:04.8 

4.7 

329 

49.99 

79.07 

Dec.  02,  1984 

03:19:06.3 

5.8 

329 

49.96 

78.86 

Dec.  16,  1984 

03:55:02.7 

6.1 

329 

mm  49.86 

78.75 

Dec.  28, 1984 

03:50:10.7 

6.0 

329 

W  49.88 

78.82 

Feb.  10, 1985 

03:27:07.6 

5.9 

329 

*  37.06 

-116.04 

Mar.  15,  1985 

16:31:00.1 

4.8 

41 

«  37.18 

-116.09 

Mar.  23,  1985 

18:30:00.0 

5.3 

41 

jS  37.09 

-116.03 

Apr.  02,  1985 

20:00:00.0 

5.7 

41 

®  37.20 

-116.21 

Apr.  06,  1985 

23:15:00.0 

4.8 

41 

„  49.92 

78.97 

Apr.  25,  1985 

00:57:06.5 

5.9 

329 

g  -22.08 

-138.90 

Apr.  30,  1985 

17:28:57.9 

4.5 

631 

W  37.25 

-116.32 

May  02,  1985 

15:20:00.0 

5.7 

41 

-21.82 

-139.05 

May  08,  1985 

20:27:58.8 

5.7 

631 

H  -22.06 

-138.86 

Jun.  03,  1985 

17:29:58.0 

5.2 

631 

1  ■  37.25 

-116.49 

Jun.  12,  1985 

15:15:00.0 

5.5 

41 

■  37.09 

-116.08 

Jun.  12,  1985 

17:30:00.0 

4.4 

41 

l  3(  49.89 

78.88 

Jun.  15, 1985 

00:57:00.7 

6.0 

329 

37.12 

-116.12 

Jun.  26,  1985 

18:03:00.0 

4.3 

41 

49.86 

78.70 

Jun.  30,  1985 

02:39:02.7 

6.0 

329 

■  65.97 

40.86 

Jul.  18,  1985 

21:14:57.4 

5.0 

724 

"  49.95 

78.83 

Jul.  20,  1985 

00:53:14.5 

5.9 

329 

49.89 

78.15 

Jul.  25,  1985 

03:11:06.5 

5.0 

329 

$  37.30 

-116.44 

Jul.  25,  1985 

14:00:00.0 

5.2 

41 

ft  37.00 

-116.04 

Aug.  17,  1985 

16:25:00.0 

4.6 

41 

37.09 

-116.00 

Sep.  27,  1985 

14:15:00.0 

4.6 

41 

■  37.21 

-116.21 

Oct.  09,  1985 

23:20:00.0 

4.2 

41 

»  37.11 

-116.12 

Oct.  16,  1985 

21:35:00.0 

4.6 

41 

-21.85 

-138.97 

Oct.  26,  1985 

16:34:58.3 

5.4 

631 

ft  -21.86 

-138.77 

Nov.  24,  1985 

16:00:58.5 

4.7 

631 

ft  -21.87 

•138.93 

Nov.  26,  1985 

17:41:58.4 

5.8 

631 

37.05 

-116.04 

Dec.  05,  1985 

15:00:00.0 

5.7 

41 

£  38  37.24 

-116.47 

Dec.  28,  1985 

19:01:00.0 

5.3 

41 

$  37.22 

•116.18 

Apr.  10,  198 

14:08:30.1 

4.9 

41 

|  37.26 

-116.44 

Apr.  22,  1986 

14:30:00.0 

5.3 

41 

^  37.10 

•116.02 

Jun.  05,  1986 

15:04:00.0 

5.3 

41 

3 

5 

Table  IV.  Nuclear  Explosions  Recorded  at  NORESS  —  (Continued) 


LAT 

LON 

DATE 

TIME 

MB 

GRN 

37.26 

-116.50 

Jun. 25, 1986 

20:27:45.1 

5.5 

41 

37.28 

-116.36 

Jul.  17,  1986 

21:00:00.0 

5.7 

41 

37.14 

-116.07 

Jul.  24,  1986 

15:05:00.0 

4.5 

41 

49.81 

78.17 

Feb.  26, 1987 

04:58:21.8 

5.4 

49.93 

78.79 

Mar.  12,  1987 

01:57:17.3 

5.4 

49.90 

78.81 

Apr.  03,  1987 

01:17:08.1 

6.2 

49.85 

78.69 

Apr.  17,  1987 

01:03:04.7 

6.0 

60.78 

56.22 

Apr.  19,  1987 

04:00:01.1 

4.5 

60.67 

56.30 

Apr.  19,  1987 

04:05:01.0 

4.4 

49.75 

77.99 

May  06,  1987 

04:02:05.6 

5.4 

i-n 


Fifnrt  1.  Top— «v«nU  wloctod  by  C«nUr  murih  «UfT;  bottom~«v«nte  •elected  by 
SAIC,  L*  JolU  group. 
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1300  N.  17th  Street,  Suite  1450 
Arlington,  VA  £££09 

ABSTRACT 

The  Center  for  Seismic  Studies  has  been  using  NORESS  data  to  test  a  number  of 
seismic  discrimination  techniques  as  DARPA  workstation  tools.  Analysis  to  date  has 
included  the  extraction  of  frequency-domain  spectral  parameters,  time-domain  autoregres¬ 
sive  parameters,  and  particle-motion  information.  The  data  set  being  used  for  testing 
includes  NORESS  recordings  of  approximately  100  events  in  Scandinavia  and  the  western 
Soviet  Union. 

Particle-motion  information  is  extracted  using  all  the  three-component  sensors  in  the 
array  by  the  covariance-matrix  averaging  technique  (Jurkevics,  1986a, b).  This  informa¬ 
tion  is  characterized  by  various  attributes  of  the  polarization  ellipse,  separately  for  dif¬ 
ferent  frequency  bands.  Power-spectral  estimates  are  obtained  from  the  autoregressive 
(AR)  modeling  of  short-period  and  high-frequency  NORESS  data,  and  the  important 
information  in  the  power  spectrum  is  condensed  into  a  few  physical  parameters.  The 
resonance  frequency  and  damping  factor  of  an  equivalent  damped  harmonic  oscillator  are 
used  to  characterize  the  spectral  content  of  the  signals  (Dysart,  1986). 

Characteristics  of  regional  phases  Pn,  Sn  and  have  been  determined  using  these 
and  other  techniques  for  blasts  and  earthquakes  in  the  area  of  interest.  For  example,  Pn 
particle  motions  tend  to  be  rectilinear  at  the  onset  of  the  phase,  and  polarization  is 
stronger  at  higher  frequencies.  For  SNR  >  2,  Pn  can  be  confidently  distinguished  from 
noise  and  L ?  on  the  basis  of  rectilinearity  and  vertical-to-horizontal  amplitude  ratio.  Sn 
particle  motion  is  also  polarized  at  the  onset  of  the  phase,  but  more  difficult  to  identify 
because  of  lower  effective  SNR.  The  main  characteristic  of  Sn  polarization  is  the  domi¬ 
nance  of  horizontal  over  vertical  motion.  The  Lf  wavetrain  is  generally  longer  and  more 
complex  than  the  other  phases,  and  the  most  consistent  Lf  motions  are  observed  around 
the  time  of  the  envelope  maximum. 

As  an  example  of  the  regional  discrimination  problem  we  analyzed  an  m^(NElS)  4.6 
event  that  occurred  on  1  August  1986,  during  the  recent  Soviet  nuclear  testing  morato¬ 
rium.  The  event  was  located  at  73.02  *  N,  56.50  *  E,  south  of  Matochkin  Shar  strait  on  the 
eastern  shore  of  Novaya  Zemlya,  U.S.S.R.  Of  34  events  that  occurred  within  200  km  of 
this  event  during  a  20-year  period  prior  to  1986,  all  except  four  appear  to  have  been 
underground  nuclear  tests  and  the  others  were  probably  aftershocks  of  one  of  the  explo¬ 
sions.  Analysis  of  the  1  August  event  using  a  number  of  discrimination  techniques  (polar¬ 
ization,  vespagram,  spectral,  cepstral  and  autoregressive  analyses)  failed  to  eliminate  the 
possibility  that  it  could  be  a  small  underground  nuclear  test.  A  search  for  depth  phases 
was  inconclusive.  GRF  and  NORESS  signals  for  this  event  had  spectral  characteristics 
that  were  not  unlike  an  4.6  presumed  underground  nuclear  explosion  in  1977,  about  80 
km  to  the  northwest.  This  illustrates  the  difficulty  of  identifying  a  well-recorded  seismic 
event,  with  magnitude  equivalent  to  a  tamped  underground  nuclear  explosion  of  a  few 
kilotons,  located  near  a  well-studied  test  site,  in  a  stable  geologic  region.  Far  greater  dif¬ 
ficulties  would  be  encountered  in  attempting  to  use  a  limited  number  of  seismic  stations 
to  detect,  locate  and  identify  much  smaller  earthquakes,  mine  blasts  and  possible  decou¬ 
pled  nuclear  explosions  in  all  parts  of  the  U.S.S.R. 
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ABSTRACT 

Assessments  of  seismic  station  networks  to  monitor  underground  nuclear  explosions 
are  usually  limited  to  estimates  of  given  station  networks.  That  is  to  say,  the  geographi¬ 
cal  coordinates  and  noise  characteristics  of  the  seismological  stations  are  given,  and  the 
capability  of  the  network  to  detect  a  seismic  event  is  then  calculated.  From  the  point  of 
view  of  monitoring,  the  reverse  problem  is  also  relevant.  That  is  to  say,  determine  the 
station  network  that  satisfies  given  desirable  performance  criteria.  This  problem  is  sel¬ 
dom  addressed,  partly  because  of  its  complexity  and  of  the  fact  that  widely  accepted  per¬ 
formance  criteria  are  not  available.  In  this  paper  we  attempt  to  get  some  insight  into  this 
problem  by  presenting  results  from  calculations  of  the  capability  of  hypothetical  networks 
with  regard  to  event  detection,  location,  and  depth  estimation. 

Calculations  for  hypothetical  networks  with  equal  station  spacing  on  a  sphere  sug¬ 
gest  that  about  10  stations  or  more  provides  fairly  uniform  detection  thresholds  over  the 
sphere,  although  there  are  significant  variations  in  areas  close  to  where  stations  are 
located.  The  calculations  also  show  that  the  average  detection  threshold  is  reduced  some¬ 
what  faster  than  the  logarithm  of  the  square  root  of  the  number  of  stations  (  log10\/n  ) 
and  appears  to  be  asymptotic  to  log1Q\/n  for  a  large  number  of  stations.  Detection 
thresholds  of  stations  located  on  shields  and  platforms  are  frequently  observed  to  be  lower 
than  those  of  stations  in  regions  of  a  different  tectonic  nature.  It  appears,  however, 
impossible  to  achieve  a  reasonably  uniform  global  coverage  from  networks  based  entirely 
on  stations  in  shields  and  platforms  areas.  Comparisons  are  also  made  for  hypothetical 
networks  employing  either  single  sensors  or  arrays.  They  suggest  that  a  network  employ¬ 
ing  small  noise  reducing  arrays  would  give  a  slightly  lower  detection  threshold  than  that 
of  a  larger  network  with  the  same  total  number  of  single  sensors. 

Location  errors  are  calculated  for  networks  with  stations  that  have  equal  distance  to 
the  epicenter  and  are  uniformly  distributed  around  the  epicenter.  The  calculated  exam¬ 
ples  suggest  the  a  90  degree  coverage  of  the  source-to-station  azimuth  gives  location  errors 
that  are  only  marginally  larger  than  those  for  networks  with  complete  azimuthal  cover¬ 
age. 

The  standard  deviation  of  the  estimated  focal  depth  on  the  basis  of  first  arrival 
times  is  used  as  a  measure  of  the  depth  determination  capability.  This  standard  devia¬ 
tion  is  entirely  determined  by  the  distance  range  covered  by  detecting  stations.  A  cover¬ 
age  limited  to  the  teleseismic  range  30-80  degrees  gives  a  comparatively  large  standard 
deviation,  whereas  a  significant  improvement  is  obtained  if  observations  at  15  degrees  or 
closer  can  also  be  utilized. 

The  results  of  the  calculations  have  been  drawn  upon  in  developing  network  concepts 
for  the  Group  of  Scientific  Experts  at  the  UN  Conference  on  Disarmament  in  Geneva. 
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Polarization  Analysis  Using  an  Array  of  Three-Component  Sensors 


ANDY  JURKEV1CS  (Center  for  Seismic  Studies,  1300  N.  17th  Street, 

Suite  1450,  Arlington,  VA  22209) 

A  processing  method  is  presented  which  computes  ground-motion 
polarization  from  an  array  of  three-component  seismometers.  The 
advantage  gained  by  using  an  array  is  similar  to  the  signal  improve¬ 
ment  realized  in  array  beamforming,  although  the  philosophy  is 
somewhat  different.  The  method  works  best  when  seismic  noise  and 
local  earth  scattering  effects  are  uncorrelated  between  array  sensors. 
In  this  case  the  variance  of  computed  particle  motions  goes  as  1/N, 
where  N  is  the  number  of  seismometers  in  the  array. 

The  analysis  is  carried  out  in  the  time  domain  and  is  very  fast.  An 
assumption  is  made  that  each  frequency  component  has  stationary 
polarization  over  several  cycles  duration.  The  data  traces  are 
bandpass  filtered  and  time  windowed  into  short  segments.  The  3x3 
covariance  matrix  is  computed  at  each  sensor,  and  the  covariance 
matrices  for  all  sensors  are  averaged  together.  The  polarization  ellip¬ 
soid  is  then  computed  using  an  eigenvalue  decomposition  and  various 
particle-motion  attributes  are  extracted  from  the  polarization  ellipse. 
The  frequency  and  time  resolutions  are  controlled  by  the  frequency 
bandwidths  and  time-segment  lengths.  The  motivation  for  the  covari¬ 
ance  matrix  averaging  comes  from  the  effectiveness  of  periodogram  or 
autocorrelation  ensemble  averaging  in  power  spectral  estimation.  The 
polarization  analysis  requires  the  velocities  of  coherent  wavefronts 
across  the  array  in  order  to  time-align  the  data  segments.  However, 
the  wavefront  velocities  are  not  required  to  nearly  the  same  accuracy 
as  needed  for  beamforming.  Beamforming  requires  a  time  alignment 
between  sensors  of  better  than  1/4  cycle  for  constructive  addition  of 
phases.  The  covariance  matrix  averaging  requires  only  that 
corresponding  time  segments  between  sensors  be  approximately 
aligned. 

Examples  are  presented  using  synthetic  data  and  seismograms 
recorded  at  NORESS,  the  small-aperture  seismic  array  in  Norway. 
The  nature  of  scattered  P  coda  is  examined,  source  aiimuths 
estimated,  and  different  regional  seismic  phases  are  distinguished  on 
the  basis  of  their  particle-motion  signatures.  This  research  is  part  of 
a  DARPA-supported  effort  to  develop  seismic  processing  capabilities 
for  Sun  workstations  at  the  Center  for  Seismic  Studies. 
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Capabilities  and  Limitations  of  the  GSETT 
Global  Seismic  Station  Network 


HANS  ISRAELSSON  (Center  for  Seismic  Studies.  1300  N  17th 
St..  Suite  1450.  Arlington.  VA  22209-3871) 

In  the  fall  of  1984  the  Ad  Hoc  Group  of  Scientific  Experts 
(GSE)  at  the  UN  Conference  of  Disarmament  carried  out  a 
world  wide  Technical  Test  (GSETT)  in  order  to  study  some 
aspects  of  an  International  Co-operative  System  for  monitoring 
a  nuclear  test  ban.  More  than  70  seismological  stations  contri¬ 
buted  measurements  over  a  period  of  two  months.  Data  from 
the  stations  were  collected  and  compiled  at  the  Center  for 
Seismic  Studies  in  Arlington,  Virginia,  which  acted  as  one  of 
three  International  Data  Centers  during  this  test.  The  col¬ 
lected  data  constitute  a  unique  set  of  measurements  in  that  it 
included  an  extended  set  of  parameters  not  ordinarily  reported 
and  many  stations  took  unusual  care  in  the  analysis  of  data. 

The  data  base  at  the  Center  in  Arlington  has  been  used  to 
examine  some  seismological  aspects  of  the  GSETT  station  net¬ 
work.  It  is  concluded  that  there  is  a  significant  variation  in  sig¬ 
nal  detection  capability  among  the  stations  (mt~4— 6),  and 
they  also  have  a  highly  non-uniform  geographical  distribution. 
This  results  in  a  network  detection  capability  which  is  strongly 
regionally  dependent  (m4« 3.C-4.6).  The  capability  of  the 
GSETT  network  varies  not  only  between  the  Northern  and  the 
Southern  Hemispheres  but  there  is  also  significant  lateral  varia¬ 
tion  within  both  Hemispheres.  The  non-uniform  geographical 
►  station  distribution  and  variation  in  station  detection  capabil¬ 

ity  also  contribute  to  the  high  percentage  of  unassociated  signal 
detections  (about  40%).  There  is  furthermore  a  significant  gap 
between  the  network  event  detection  thresholds  and  thresholds 
for  carrying  out  other  major  functions.  For  example,  the  loca¬ 
tion  accuracy  of  the  GSETT  network  decreases  significantly  at 
mk~5.  The  capability  of  the  GSETT  network  to  determine 
surface  wave  magnitudes  and  compile  other  identification 
parameters  is  also  significantly  lower  than  that  for  ev-nt  detec¬ 
tion. 

Computer  experiments  with  a  network  simulation  program 
(SNAP/DI  have  been  carried  out  to  study  how  the  the  GSETT 
network  should  be  supplemented  and  upgraded  so  that  it  will 
perform  in  a  more  uniform  and  balanced  manner. 


Autoregressive  Analysis  of  Regional  Earthquakes  and  Explosions 
PAUL  S.  DYSART  (Center  for  Seismic  Studies,  Rossiyn,  VA  22209) 

The  prediction  coefficients  in  an  autoregressive  (AR)  model  contain  ampli¬ 
tude  and  frequency  information  which  is  easily  and  quickly  obtained 
without  Fourier  transforming  the  time  series.  By  relating  the  prediction 
coefficients  of  the  second-order  AR  model  to  a  damped  harmonic  oscillator, 
this  representation  of  seismic  signals  yields  a  set  of  spectral  parameters 
which  are  very  useful  in  problems  of  regional  event  discrimination,  detec¬ 
tion,  and  phase  identification.  Variations  in  the  AR  parameters  with  time 
reflect  changes  in  the  spectral  character  of  seismic  signals  without  the 
necessity  of  decomposing  the  seismogram  into  separate  passbands.  When 
viewed  in  this  manner  the  AR  parameters  provide  a  promising  means  of 
detecting  and  identifying  Pn,  Lg,  and  other  regional  phases.  An  AR  algo¬ 
rithm  has  been  developed  as  part  of  the  design  of  a  DARPA  seismic  works¬ 
tation.  In  most  cases  the  speed  of  computation  allows  an  automated  spec¬ 
tral  analysis  in  real  time. 

In  this  study,  a  resonance  frequency,  damping  factor,  power,  and  final 
prediction  error  (FPE)  derived  from  the  second-order  AR  model  coeficients 
are  used  to  characterize  regional  earthquakes  and  mine  blasts  recorded  by 
the  NORESS  array  in  southern  Norway.  These  parameters  are  computed 
as  functions  of  time  with  frequency  and  time  resolutions  controlled  by  the 
length  and  overlap  of  a  sliding  window.  Results  indicate  that  the  Pn  reso¬ 
nance  frequency,  which  estimates  the  dominant  frequency  of  the  spectrum, 
is  higher  for  the  explosions,  as  would  be  expected  from  events  with  shorter 
source  duration.  The  damping  factor  or  bandwidth  is  greater  for  explo¬ 
sions,  as  is  the  FPE,  which  is  shown  to  be  a  good  estimate  of  the  signal 
entropy.  As  discriminants,  both  the  resonance  frequency  and  the  FPE  for 
Pn  waves  separate  the  explosions  from  the  earthquakes,  although  there  is 
an  overlapping  of  the  two  populations.  It  appears  that  some  correction  for 
attenuation  with  distance  is  necessary  for  the  more  distant  events.  The 
NORESS  results  are  consistent  with  the  AR  analysis  of  Eurasian  events 
recorded  at  NORSAR  (Tjostheim  et  al.,  1975).  A  comparison  of  the  two 
studies  implies  that  the  more  extensive  set  of  Eurasian  events  also  exhibits 
higher  frequencies  and  wider  bandwidths  for  underground  nuclear  explo¬ 
sions  than  for  earthquakes. 


The  Character  of  Regional  Swwiwfrum  Recorded  by  the  NORESS  A rrej 
ImptkitioM  for  E»wt  MmlifiutiM  and  Automata*  Procettmg 

JAY  J  PULL! 

PAUL  S  DYSART  (both  at  Center  for  Seismic  Studies  Rottiyn  VA  22709) 

NORESS  it  a  25  element  small  apart  ura  mwik  array  toe  at  a*  m  southern 
Norway  Tha  array  currently  utat  an  automata*  procettmg  tyttam  for  tha 
detection  art*  location  of  rational  avantt  an*  na war  tyttami  ara  un*ar  develop 
mant  which  will  incorporate  artificial  intelligence  (Al)  technologies  Tha  goal  of 
our  r  eta  arch  it  to  invattigata  tha  waveform  characteristics  of  avantt  recor*e*  at 
NORESS  an*  datarmino  whether  or  not  theta  characteristics  can  ba  uto*  for 
avant  identification  by  an  Al  tyttam 

We  have  analyze*  approximately  SO  earthquakes  an*  chemical  upteueei 
which  occurred  within  1500  km  of  NORESS  Dittinct  an*  repeatable  waveform 
cha/ act  eristics  can  ba  observed  for  groupt  of  avantt  Tor  an  ample  Lg  wave* 
generate*  by  chemical  explosions  ara  ban*  km  it  a*  to  latt  than  I  Hz  whereat 
Lg  wavat  generate*  by  continental  earthquakes  often  extan*  to  higher  frequen 
ciat.  Tha  largatt  waveform  variationi  by  far  ara  for  Pn  wavat  Pn  wavat  gen 
crate*  by  many  explotione  in  tha  wettem  Soviet  Union  arrive  at  NORESS  at  a 
wave  packet  with  a  tharply  peake*  tpectrum  Theta  tpaclra  alto  thow  scallop 
ing  at  harmonic  fraquencict  which  may  ba  indicative  of  delayed  or  multiple 
shoti.  Thit  it  alto  evident  in  tha  cepetraf  analytic  of  Pn  Pn  wavat  for  earth 
quakes  at  timilar  distance!  arc  recorded  at  dittinct  multiple  arrivalt  Many 
explosions  in  Scandinavia  generate  Pn  wavat  with  peak  fraquencict  in  tha  t  to  • 
Hz  range,  whereas  the  Pn  wavat  for  earthquakes  of  timilar  magnitudes  ara  of 
lower  frequency. 

Our  experience  indicates  that  no  tingle  rule  or  tat  of  rules  can  be  universally 
applied  at  NORESS  to  identify  source  types  Such  an  approach  would  require 
that  path  corrections  ba  known  to  a  high  degree  of  accuracy  However  tha 
repeatability  of  waveform  characteristics  observed  with  tha  NORESS  array  tug 
gests  that  a  large  amount  of  information  may  ba  accumulated  in  a  short  period 
of  time,  and  event  identification  may  be  greatly  aided  by  applying  rare- freer* 
rather  than  rule- bated  reasoning. 


S*i»mic  Discrimination  and  the  1  August  1986  Novaya  Zemlya  Event 

A  S  RVALL,  A  JURKEVICS,  P  S  DYSART  and  J  J  PULLI  (all  at  Ctr. 

for  Seismic  Studies,  1300  N.  17th  St.,  Suite  1450,  Arlington,  VA 

22209) 

As  an  example  of  the  difficulty  of  the  discrimination  problem  in 
nuclear  test  ban  treaty  verification,  we  analysed  an  4.6  event  that 
occurred  on  1  August  1996,  during  the  recent  Soviet  nuclear  testing 
moratorium.  The  event  was  located  at  73.02  *  N,  56.50  *  E,  south  of 
Matochkin  Shar  strait  on  the  eastern  shore  of  Novaya  Zemlya,  USSR. 
Of  34  events  that  occurred  within  200  km  of  this  event  during  a  20- 
year  period  prior  to  1966,  all  except  four  appear  to  have  been  under¬ 
ground  nuclear  tests  and  the  others  were  probably  aftershocks  of  one 
of  the  explosions.  Analysis  of  the  1  August  event  using  a  number  of 
discrimination  techniques  being  tested  as  DARPA  workstation  tools 
at  the  Center  for  Seismic  Studies  (polarisation,  vespagram,  spectral, 
cepstral  and  autoregressive  analyses)  failed  to  eliminate  the  possibil¬ 
ity  that  it  could  be  an  underground  nuclear  teat.  A  search  for  depth 
phases  was  inconclusive.  Signals  recorded  for  this  event  by  the  GRF 
and  NORESS  arrays  had  spectral  characteristics  that  were  not  unlike 
an  m,  4.6  presumed  underground  nuclear  explosion  in  1977,  about  80 
km  to  the  northwest. 

These  results  illustrate  the  difficulty  of  identifying  a  well-recorded 
seismic  event,  with  magnitude  equivalent  to  a  tamped  underground 
nuclear  explosion  of  a  few  kilotons,  located  near  a  well-studied  test 
site,  in  a  stable  geologic  region.  Far  greater  difficulties  would  be 
encountered  in  attempting  to  use  a  limited  number  of  seismic  stations 
to  detect,  locate  and  identify  much  smaller  earthquakes,  mine  blasts 
and  possible  decoupled  nuclear  explosions  in  all  parts  of  the  USSR. 


